
Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 


1. Thesis and Dissertation Collection, all items 


2002-09 

Remotely sensed density measurements of 
volcanic sulfur dioxide plumes using a 
spectral long wave infrared imager 

Mares, Aimee Gail 

Monterey, California. Naval Postgraduate School 


http://hdl.handle.net/10945/4879 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 


Downloaded from NPS Archive: Calhoun 



DUDLEY 

KNOX 

LIBRARY 


http ://w w w. nps.edu/lEbrary 


Calhoun is the Naval Postgraduate Schools public access digital repository for 
research materials and institutional publications created by tire NPS community, 
Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 
appointed — and published — scholar^ author. 

Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 


NAVAL POSTGRADUATE SCHOOL 
Monterey, California 



THESIS 


REMOTELY SENSED DENSITY MEASUREMENTS OF 

VOLCANIC SULFUR DIOXIDE PLUMES USING A 

SPECTRAL LONG WAVE INFRARED IMAGER 

by 


Aimee Gail Mares 


September 2002 


Thesis Advisor: 

Richard C. Olsen 

Second Reader: 

Paul G. Lucey 


Approved for public release; distribution unlimited 






THIS PAGE INTENTIONALLY LEFT BLANK 



REPORT DOCUMENTATION PAGE 


Form Approved OMB No. 0704-0188 
Public reporting burden for this collection of information is estimated to average 1 hour per response, including 
the time for reviewing instruction, searching existing data sources, gathering and maintaining the data needed, and 
completing and reviewing the collection of information. Send comments regarding this burden estimate or any 
other aspect of this collection of information, including suggestions for reducing this burden, to Washington 
headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson Davis Highway, Suite 
1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project 
(0704-0188) Washington DC 20503. 

I. AGENCY USE ONLY (Leave blank) 2. REPORT DATE 3. REPORT TYPE AND DATES COVERED 

_ September 2002 _ Master’s Thesis _ 

4. TITLE AND SUBTITLE: Remotely Sensed Density Measurements of Volcanic 5. FUNDING NUMBERS 
Sulfur Dioxide Plumes Using a Spectral Long Wave Infrared Imager _ 

6. AUTHOR(S) Aimee Gail Mares _ 

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 8. PERFORMING 

Naval Postgraduate School ORGANIZATION REPORT 

Monterey, CA 93943-5000 NUMBER 

9. SPONSORING /MONITORING AGENCY NAME(S) AND ADDRESS(ES) 10. SPONSORING/MONITORING 
N/A AGENCY REPORT NUMBER 

II. SUPPLEMENTARY NOTES The views expressed in this thesis are those of the author and do not reflect the official 

policy or position of the Department of Defense or the U.S. Government. _ 

12a. DISTRIBUTION / AVAILABILITY STATEMENT 12b. DISTRIBUTION CODE 

Approved for public release; distribution unlimited _ 

13. ABSTRACT (maximum 200 words) 

This thesis examines the process of detecting and quantifying volcanic S0 2 plumes using the Airborne Hyperspectral 
Infrared Imager (AHI) developed by the University of Hawaii. AHI was flown over Pu’u’O’o Vent of Kilauea Volcano in 
Hawaii to collect data on S0 2 plumes. In conjunction with these observations, data were taken with the Hawaii Volcano 
Observatory’s Correlation Spectrometer (COSPEC) and University of Hawaii’s FLYSPEC. These are ultraviolet remote 
sensors with a successful history of monitoring volcanic S0 2 plumes at ~0.3 p. 

AHI is a LWIR pushbroom imager sensitive to the 7.5 - 11.5 p region. Spectral analysis and mapping tools were 
used to identify and classify the S0 2 plume in both radiance and emissive space. MODTRAN was used to model the radiance 
observed by the sensor as it looked to the ground through an S0 2 plume. A spectral library of radiance profiles with varying 
ground surface temperatures and S0 2 concentrations was developed, and the AHI data fitted to the varying model profiles. 
Reasonable values of S0 2 emission were obtained, though the values directly over the vent obtained by AHI were much higher 
than those obtained by the UV sensors some distance away. 

14. SUBJECT TERMS hypersepctral, airborne remote sensing, Long Wave Infrared 15. NUMBER OF 
(LWIR), Ultraviolet (UV), Sulfur Dioxide, gas detection, volcanic plumes, Airborne Hyperspectral PAGES 
Infrared Imager (AHI), FLYSPEC 146 

16. PRICE CODE 

17. SECURITY 18. SECURITY 19. SECURITY 20. LIMITATION 

CLASSIFICATION OF CLASSIFICATION OF THIS CLASSIFICATION OF OF ABSTRACT 

REPORT PAGE ABSTRACT 

Unclassified Unclassified Unclassified UL 


NSN 7540-01-280-5500 


Standard Form 298 (Rev. 2-89) 
Prescribed by ANSI Std. 239-18 


1 




























THIS PAGE INTENTIONALLY LEFT BLANK 


11 



Approved for public release; distribution unlimited 


REMOTELY SENSED DENSITY MEASUREMENTS OF VOLCANIC SULFUR 
DIOXIDE PLUMES USING A SPECTRAL LONG WAVE INFRARED IMAGER 


Aimee G. Mares 

Captain, United States Marine Corps 
B.S., United States Naval Academy, 1995 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE IN SPACE SYSTEMS OPERATIONS 


from the 


NAVAL POSTGRADUATE SCHOOL 
September 2002 


Author: Aimee Gail Mares 


Approved by: Richard C. Olsen 

Thesis Advisor 


Paul G. Lucey 
Second Reader 


Rudy Panholzer 

Chairman, Space Systems Academic Group 



THIS PAGE INTENTIONALLY LEFT BLANK 


IV 



ABSTRACT 


This thesis examines the process of detecting and quantifying volcanic SO 2 
plumes using the Airborne Hyperspectral Infrared Imager (AHI) developed by the 
University of Hawaii. AHI was flown over Pu’u’O’o Vent of Kilauea Volcano in Hawaii 
to collect data on SO 2 plumes. In conjunction with these observations, data were taken 
with the Hawaii Volcano Observatory’s Correlation Spectrometer (COSPEC) and 
University of Hawaii’s FLYSPEC. These are ultraviolet remote sensors with a successful 
history of monitoring volcanic SO 2 plumes at ~0.3 p. 

AHI is a LWIR pushbroom imager sensitive to the 7.5 - 11.5 p region. Spectral 
analysis and mapping tools were used to identify and classify the SO 2 plume in both 
radiance and emissive space. MODTRAN was used to model the radiance observed by 
the sensor as it looked to the ground through an SO 2 plume. A spectral library of 
radiance profiles with varying ground surface temperatures and SO 2 concentrations was 
developed, and the AHI data fitted to the varying model profiles. Reasonable values of 
SO 2 emission were obtained, though the values directly over the vent obtained by AHI 
were much higher than those obtained by the UV sensors some distance away. 
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I. INTRODUCTION 


The modem science of volcanology originated in Hawaii with the establishment 
of the Hawaiian Volcano Observatory (HVO) in 1912. Prior to the advent of 
volcanology, volcanoes were generally only studied after a major eruption. In the early 
1900’s, Thomas A. Jaggar, Jr., a geologist at the Massachusetts Institute of Technology 
(MIT) started promoting a continuous study of volcanoes. He believed that scientists 
would gain a better understanding of volcanoes if studies were conducted before, during, 
and after eruptions. His beliefs compelled the residents of Hawaii to establish the 
Hawaiian Volcano Research Association (HVRA). With assistance from HVRA and the 
Whitney Fund of MIT, Jagger established the HVO. Today, the HVO is operated by the 
U.S. Geological Survey (USGS). (Tilling, R., Heliker, C., and Wright, T., “Eruptions of 
Hawaiian Volcanoes”) 

The mission of HVO is to monitor Hawaii’s Mauna Loa and Kilauea Volcanoes. 
“The term volcano monitoring refers to the observations and measurements scientists 
make to document changes in the state of the volcano during and between eruptions.” 
(Watson, “Eruptions of Hawaiian Volcanoes”) Through this continuous and close 
monitoring, the scientists can help protect the people of Hawaii and their property by 
issuing timely warnings of hazardous activity. Lava flows and explosive eruptions are 
the hazards most often associated with volcanoes, but additional hazards exist such as 
volcanic smog, earthquakes, and tsunamis. (Heliker, C., Stauffer, P, and Hendley, J., 
“Living on Active Volcanoes”) 

Hawaii’s Kilauea Volcano is unique in its long-term (1983 - present), nearly 
continuous eruptive activity and has provided scientists an ideal site to develop a greater 
understanding of many volcanic life processes. Another unique aspect of Kilauea 
Volcano is its approachable lava flows. Kilauea’s relatively gentle nature and convenient 
accessibility can be deceiving, as these characteristics do not eliminate any dangers 
associated with studying or viewing the volcano. 

Scientists use a wide variety of instruments, sensors, methods, and techniques to 
monitor and study volcanic activity. In order to conduct safe and continuous monitoring, 
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scientists are researching sensors and systems that will allow them to remotely collect 
data for their studies. Remote collection not only ensures safety, but it often offers a 
synoptic and repetitive view of the volcanoes that may not always be possible through 
direct collection. Mouginis-Mark et al.’s Remote Sensing of Active Volcanism is a 
compilation of many scientific papers addressing a wide range of airborne and 
spacebome remote sensing instruments that have been applied to aid in the study of 
volcanoes throughout the world. A summary of these instruments and the volcanic 
processes studied is found in Table 1. 


Instrument 

Total Ozone Mapping Spectrometer (TOMS) 

Sensor Platform 
Spaceborne 

Wavelength Range 
.300 - .340 microns 

Volcanic Process Studied 
Stratospheric sulfur dioxide 

Advanced Very High Resolution Radiometer (AVHRR) 

Spaceborne 

.58 -1.1 microns 

10.3 -12.5 microns 

Daytime identification of eruption p 
Identification of ash clouds 

Land Satellite (Landsat) 

Spaceborne 

.8 - 2.2 microns 

10.4 -12.4 microns 

Determination of temperatures and 
Detection of low temperature them 

Airborne Topographic Laser Altimeter System (ATLAS) 

Airborne 

1.06 microns 

Topographic mapping via an active 

Fourier Transform Infrared Spectroscopy (FTIR) 

Ground 

3 - 5 microns 

8-12.5 microns 

Monitoring volcanic gas emissions 
Monitoring volcanic gas emissions 

Geostationary Operational Environmental Satellite (GOES) 

Spaceborne 

3.78 - 4.03 microns 

Real-time identification of hot spot; 

Thermal Inrared Multispectral Scanner (TIMS) 

Airborne 

8-12 microns 

Retrieval of sulfur dioxide emission 

High Resolution Sounder (HIRS/2) 

Spaceborne 

8.16 -13.97 microns 

Properties of silicate ash in clouds 

Earth Resource Satellite (ERS-2), Shuttle Imaging Radar (SIR-C), and 
Topographic Synthetic Aperture Radar (TOPSAR) 

Spaceborne 

5.6 - 24 cm 

Topographic mapping using radar i 

ERS-2, Japanese Earth Resource Satellite (JERS-1) 

Spaceborne 

5.6 - 24 cm 

Deformation of volcanoes using ra< 


Table 1. Remote Sensing Instruments Used In Volcanic Monitoring 
(after Mouginis-Mark, et al., Appendix 3., p. 269) 


This paper focuses on monitoring the gas emission process of Kilauea Volcano. 
During periods of sustained eruption, Kilauea emits about 2,000 tons of sulfur dioxide 
gas (SO 2 ) each day. There are many reasons for monitoring the amount of SO 2 emitted 
from Kilauea. SO 2 is a poisonous gas that can irritate human respiratory systems. It also 
reacts with oxygen and water in the atmosphere to produce volcanic smog (vog) and acid 
rain. (Sutton, et al., “Volcanic Air Pollution”) In addition to better understanding the 
affects it may have on the health and welfare of the people of Hawaii, the SO 2 emission 
rate is also monitored with the belief that it can also aid in predicting periods of 
increasing or decreasing activity within the volcano. Sutton et al. investigate this topic in 
their study of the “implications for eruptive processes as indicated by sulfur dioxide 
emissions”. 
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The Correlation Spectrometer (COSPEC) is the accepted and standard instrument 
used by HYO and volcano observatories throughout the world in monitoring volcanic 
SO 2 emission rates. Chapter II describes the sensor, the data collection and analysis 
techniques, along with a summary of data collected for Kilauea from 1979 through 1997. 
FLYSPEC, a miniature correlation spectrometer developed by the University of Hawaii’s 
Hawaii Institute of Geophysics and Planetology (HIGP) is also described in Chapter II. 
FLYSPEC’s original design applies new and current technologies to detecting SO 2 and 
producing emission rates in a manner similar to that of the older technologies of 
COSPEC. The two instruments have collected multiple data sets simultaneously. Results 
from these data collections are also discussed in Chapter II. COSPEC and FLYSPEC, 
operating in the ultraviolet (UV) range of the electromagnetic spectrum and most often 
employed from the ground, do have some limitations. This paper looks to imaging 
airborne sensors in the infrared (IR) region of the electromagnetic spectrum that may be 
used to augment COSPEC or FLYSPEC and fill any void that may be left by these 
sensors. 

NASA’s Thermal Infrared Multispectral Scanner (TIMS) is a long wave infrared 
(LWIR) imaging sensor that has collected data on Kilauea Volcano and demonstrated the 
ability detect the SO 2 gas and produce an SO 2 emission rate. The TIMS sensor, its data 
collection and analysis techniques, and a summary of the data results are discussed in 
Chapter III. Chapter III also briefly discusses the application of TIMS analysis 
techniques to data collected by the Advanced Spacebome Thermal Emission Reflection 
Radiometer (ASTER). ASTER is another multispectral LWIR imager employed on the 
Terra Satellite as part of NASA’s Earth Observing System (EOS). A third LWIR sensor, 
the Airborne Hyperspectral Imager, is also introduced in Chapter III. The University of 
Hawaii’s HIGP originally designed and built (AHI) for detection of buried mines, but 
many additional applications are being explored. The sensor, its data collection 
techniques, and its current applications are disused in Chapter III. 

Chapter IV explores the ability of AHI to detect SO 2 and produce an emission 

rate. AHI, COSPEC, and FLYSPEC collected data on the SO 2 plume of Kilauea on April 

18, 2002. The site of the data collection, analysis techniques applied to AHI’s data, and 

the results are discussed in Chapter IV. Chapter V discusses these results and provides a 

3 



conclusion on the potential use of AHI for SO 2 detection and producing SO 2 emission 
rates. 
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II. ULTRAVIOLET DETECTION OF S0 2 

A. ULTRAVIOLET SESNORS 

The Correlation Spectrometer (COSPEC) is the proven and trusted sensor widely 
used within the volcanological community for detecting and measuring S0 2 
concentrations within volcanic plumes and determining S0 2 emission rates. COSPEC 
has been used for over 30 years, yet there has been very little advancement in its 
technology or data analysis methods. FLYSPEC is currently a developmental sensor 
using new spectrometer and computing technologies that may prove capable of detecting 
and measuring S0 2 concentrations and determining emission rates with the same 
accuracy and reliability as COSPEC, but with a much smaller system that has more 
processing capability and is less expensive to acquire, develop, and maintain. FLYSPEC 
will hopefully prove that in some cases better, faster, cheaper, and smaller is actually 
possible. 

Both sensors operate in the ultraviolet (UV) region of the electromagnetic 
spectrum and use similar concepts for detection, measurement, and data processing. The 
following paragraphs provide an overview of the system and sensor designs, the data 
collection methods, and basic calculation methods used by COSPEC and FLYSPEC. 

1. COSPEC 

Barringer Research, Toronto, Canada, developed COSPEC primarily for 
environmental monitoring of S0 2 in the 1960s. Researchers first used COSPEC to 
monitor the S0 2 flux from volcanoes in April 1971, at Mt. Mihara in Japan. The use of 
COSPEC for volcanic S0 2 flux measurements spread throughout the world providing a 
collection of data that has been used to estimate a worldwide S0 2 flux from volcanoes, 
delineate the shapes and concentrations within volcanic plumes, and to determine 
temporal variations of S0 2 flux from volcanoes. (Tazieff and Sabroux, 1983, p. 425) 
a. System Design 

The COSPEC system can be described in two segments: a data collection 
segment and a data analysis segment. These segments are sequential, as COSPEC does 
not have the ability to process and analyze the data as it is collected. 
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(1) Data collection segment. The data collection segment is 
normally composed of three main components: the UV sensor, a portable computer or 
data logger, and an analogue paper chart recorder as illustrated in Figure 1. Solar UV 



Figure 1. COSPEC System Components of Data Collection Segment 


radiation passes through the UV sensor, which is connected to both the data logger and 
paper chart recorder. The paper chart recorder displays analog voltage recordings, which 
correspond to a real time SO 2 absorption profile for the area covered by the sensor. The 
data logger records these voltage readings in a digital format. Figure 2 summarizes the 
data collection segment including inputs and outputs of the system. 



UV Sensor 


\ 


Portable Computer 
And / or 
Data logger 


Analogue 
Paper chart 
Recorder 


Output = Digital and analog voltage 
Recordings which correspond to 
S02 absorption profiles / spectrum 



Figure 2. COSPEC System Inputs, Components, and Outputs Used for 

Data Collection Segment 

(2) Data analysis segment. The digital voltage recordings, 

collected in the field, are transferred from the data logger to a desktop computer. The 

data processing produces an average SO 2 concentration path length with units of parts per 

million-meter (ppm-m). This average concentration path length is then combined with 

wind data (speed and direction), location data, and time to produce an average SO 2 

emission rate for the day with units of tons per day (t/d). The data analysis segment and 
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the information it produces are illustrated in Figure 3. Details of this data processing and 
calculations are given in Section d. 


U Input = Digital 
voltage recordings 

Figure 3. COSPEC System Inputs, Components, and Outputs of the 

Data Analysis Segment 

b. Sensor Design 

Resonance Inc. is the company currently manufacturing COSPEC sensors. 
Several physical and performance specifications for the sensor are listed in Table 2. 
COSPEC’s sensor design is normally described in three sections: front optics, mid 
optics, and rear optics. (Barringer, 1976, p. 7) These sections, the components that 
comprise them, and the general flow of radiation through the sensor are depicted in 
Figure 4. 

(1) Front Optics. Solar radiation is collected by the Cassegrain 
style telescope, which has a field of view of 23 milliradians by 7 milliradians (1.3178 
degrees by .401 degrees). This collected radiation passes through an entrance slit, which 
reduces the amount of radiation that will enter the mid optics section. 

(2) Mid Optics. The radiation is directed and focused onto a 
diffraction grating, by a series of mirrors in the mid section of the sensor. The diffraction 
grating separates the radiation into individual wavelengths, and this dispersed radiation is 
directed toward the rear optics section. 

(3) Rear Optics. In the rear section of the sensor, the radiation 
is focused on a correlator disc, which has arrays of circular slits etched in it correlating to 
positive and negative (peaks and troughs) SO 2 absorption bands. A photomultiplier tube 
monitors the radiation modulated by these slits. Changes in the ratio of radiation coming 
through the slits indicate a presence or absence of SO 2 and are proportional to the 
concentration of SO 2 detected. These ratio changes are reflected in the output voltages 
processed by the electronics and displayed in a text format on the data logger and as a 
real time absorption profile on the paper chart recorder. It should be noted that the units 


Desktop Computer for 
Data Processing / Reduction 
To produce average 
Concentration length 


A I 
A| 


Permanent Wind Station 
approx. 3 meters above 
Ground to measure wind speed 


Map and accurate time source to determine 
specific location at specific times along a 
predetermined “straight line path 


Output = 

S02 emssion rate , 
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of measure of the real time absorption profile are simply volts and require post¬ 
processing to obtain SO 2 column content. 

The rear optics also includes two quartz glass calibration cells. 
Each cell contains a known concentration of SO 2 gas. One of the cells contains a 
concentration that is anticipated to be relatively higher than any actual SO 2 concentration 
measured, and the other cell contains a concentration that is anticipated to be relatively 
low in comparison to what is being measured. The calibration cells aid in determining 
sensor response and detecting any sensor drift that may be present. (Tazieff and Sabroux, 
1983, p. 427) 



Incoming solar radiation 

I 

| Cassegrain Telescope (1) | 

t _ 

Entrance Slit (2) 

Reduces entrance radiation 

t 

Mirror (3) 

Directs radiation 

_ t 

Spherical Mirror (4) 

Focuses radiation 

t 

Diffraction Grating (5) 

Separates radiation into 
individual wavelengths 

t 

| Mirror (6) | 

_t_ 

Correlator Disc (7) 

Allows passage of only two sets of wavelengths 
through slits. One set corresponds to S02 
absorption And the other to no S02 absorption. 

t 

Calibration Cells (8) 

t 

Photomultiplier Tube (9) 

,_ * _ , 


^ Front Optics 




V Mid Optics 


J 


Rear Optics 


Output to data logger and 
analog paper chart recorder 


Figure 4. COSPEC Sensor Components (After Barringer, Figure 1, p. 8) 


c. Data Collection Techniques 


COSPEC has three standard techniques of collecting data on volcanic 
plumes: airborne, ground stationary, and ground mobile. Each technique has advantages 
and disadvantages that revolve around the location, access to, surrounding terrain, and 
environment of the volcano. Typical plume size, shape, and altitude also effect decisions 
as to which collection technique may be best. This presentation will focus on the ground 
mobile technique. 


Using the ground mobile technique, COSPEC is typically mounted in a 
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vehicle as illustrated in Figure 1. The telescope is pointed out the passenger side window 
with a clear view of the sky directly above it. In an ideal scenario, the vehicle will be 
able to make several traverses under the volcano’s plume on roads that are close to the 
vent and as perpendicular to the direction of the plume drift as possible. The traverses 
are normally broken down into segments along the road, which are relatively straight 
with uniform geometry to the plume. These road segments are manually annotated on the 
absorption profile generated by the paper chart recorder as the traverses are made. 
During data analysis, data from each of these segments is processed separately. This 
allows individual geometry corrections to be applied to any road segments which may not 
have been exactly perpendicular to the plume. Figure 5 illustrates an ideal scenario for 
using the ground mobile technique. (Tazieff and Sabroux, 1983, p. 433) 

As Figure 5 illustrates, one traverse is typically broken into several 
straight-line segments. The plume will not be continuously visible for the entire traverse. 
The vehicle will be driving in and out of the plume. It is standard procedure to start 
collecting data outside of the plume under “clean air”, air that is free of SO 2 . At the 
beginning of each traverse, calibration data is collected using both the high and low 
calibration cells. This is accomplished by rotating the calibration cells into the sensor’s 
field of view. It is very important to ensure the vehicle is in “clean air” as the calibration 
data is recorded. This is fairly easy to accomplish using the real time absorption profile 
generated by the analog paper chart recorder. Once the calibration data is collected, the 
sensor is driven through the traverse. At the end of the traverse, a collection of 
calibration data is repeated. 



Figure 5. Ideal Scenario for Collecting Data Using the Ground Mobile 

Technique 
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The ground mobile technique works best when the plume is “on the 
ground”, and when there is convenient road access with appropriate geometry close to the 
vent and plume. It is logistically easier, requires fewer people, and is less costly with 
respect to both money and time when compared to other techniques. (Tazieff and 
Sabroux, 1983, p. 433) 

d. Data Analysis and Calculations 

COSPEC has been monitoring volcanic SO 2 plumes for 30 years. 
Extensive records of SO 2 emission rates have been maintained for volcanoes throughout 
the world. The upkeep of these records and establishing a historical data set is an 
important aspect in improving the understanding of volcanic processes and their effects 
on the environment. In an effort to maintain consistency over time and throughout the 
world, two standard calculations for analysis of COSPEC data have been established to 
compliment the standard data collection techniques. 

(1) Calculation of Average Concentration Path Length. The 
first calculation converts the SO 2 plume’s absorption profile values, initially in units of 
volts, to an average concentration path length in parts per million-meter (ppm-m) for one 
segment of the traverse. 

One ppm-m of SO 2 is equivalent to one cubic centimeter of SO 2 gas 
uniformly mixed in one million cubic centimeters of air that is viewed by 
COSPEC over an optical path of one meter at a pressure of 101.325 kPa 
and a temperature of 20 degrees Celsius. (Jones and Stix, 2001, p. 44) 

The following equation is used to calculate this average concentration path length per 
segment in ppm-m: 

SChseg= — * SChca, (1.1) 

Peal 

S 02 se g = average concentration path length in ppm-m 
S 02 cai = known concentration of the appropriate calibration cell in ppm-m 
Peal = peak height of the appropriate calibration cell in arbitrary units 
P seg = average segment peak height also in arbitrary units 

P seg must be calculated from the following equation in order to choose the appropriate 
values for the P ca i and S02 Ca i variables in Equation 1.1: 
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_ Segment Area 

Pseg — - (1.2) 

Segment Width 

If P seg is less than the low calibration peak height, then the low 
calibration cell values should be used for the S02 Ca i and P ca i variables in equation 1.1. If 
P seg falls between the high and low calibration heights, then an average of the two 
calibration cell values should be used for the SC^cai and P ca i variables. Finally, if P seg is 
greater than the high calibration peak height, then the high calibration cell value should 
be used for the S02 Ca i and P ca i variables. A simulated COSPEC absorption profile is 
shown in Figure 6. In this figure P seg is between the high and low calibration cell heights, 
so the average of the calibration cell values would be used for calculating Equation 1.1. 
(Jones and Stix, 2001, p. 45) 



Figure 6. Simulated COSPEC Absorption Profile Illustrating Variables 

Used in the Concentration Path Length Calculation in Equation 1.1 

(2) Calculation of Emission Rates. The second calculation 
determines an emission rate in tons per day (t/d) by taking the average concentration path 
length calculated in Equation 1.1 and combining it with wind data (speed and direction), 
segment length, and any angle corrections required to create a perpendicular geometry 
between the plume and road segment. An emission rate is calculated for each segment. 
These values are then summed to determine the total SO 2 emission rate for a given 
traverse. The SO 2 emission rate calculation for each segment is shown in Equation 1.3. 
(Jones and Stix, 2001, p. 44) 

Eseg = SOlseg * COS 0 * dl * Vwind * Cf (1.3) 

E = SO 2 emission rate per segment in tons per day (t/d) 

S 02 seg = average concentration path length per segment in ppm-m 

cos 0 = angle correction required to create perpendicular geometry between the 
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plume and road segment 

dl = length of a segment in meters determined from a map 

v W ind = average wind speed in m/s determined by wind station or other means 

Cf = conversion factor to convert from ppm-m 3 /s to metric tons per day (t/d) 
Equation 1.4 calculates the conversion factor: 

Cf = 2.8579 * .001 * 10 _6 * 86400 *273/293 = .00023008194 (1.4) 

2.8579 = density of SO 2 gas in kg/m at standard temperature (0 degrees Celsius) 
and pressure (101.325 kPa) 

.001 = conversion of kg to g 

10~ 6 = conversion of g to tons 

86400 = conversion of seconds to days 

273/279 = converts SO 2 gas density from standard temperature (0 degrees 
Celsius) to 20 degrees Celsius 

These collection techniques and analysis methods have faithfully generated trusted data 
calculations for three decades; however, this does not mean COSPEC is a perfect system. 
As with many systems, advancements in technology create opportunities for 
improvements that should not be overlooked. 

2. FLYSPEC 

Keith Horton, Assistant Researcher at the Hawaii Institute of Geophysics and 
Planetology, is currently developing a sensor referred to as FLYSPEC (a name given to 
emphasize its size in comparison to COSPEC). FLYSPEC is essentially a miniature 
COSPEC that takes advantage of new spectrometer and computing technology. The new 
technology reduces the size and cost of the system, in addition to reducing the amount of 
time required to manually analyze the data. More importantly, the new technology offers 
real time data analysis capability with spectral fitting algorithms that could prove to 
provide more accurate results. 

a. System Design 

As FLYSPEC is a developmental sensor, no formal literature is currently 
published. For comparison purposes, the FLYSPEC system can be described in the same 
two segments as COSPEC: a data collection segment and a data analysis segment. With 
FLYSPEC these segments are not strictly sequential. The data collection segment does 
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perform some data processing and produces real time data that can be read in ppm. This 
difference is explained in further detail in the following paragraphs. 


(1) Data Collection Segment. FLYSPEC’s data collection 


segment is composed of three main components: the UV sensor, a portable computer, and 
a Global Positioning System (GPS). These three components are pictured in Figure 7. 
Solar UV radiation passes through the UV sensor which is connected to the portable 



parallel cable connection option 


Figure 7. FLYSPEC System Components Used in the Data 

Collection Segment 


computer via either a USB or a parallel port cable. The portable computer displays the 
real time UV spectrum seen by the sensor along with a real time SO 2 absorption profile 
for the area covered by the sensor. To accompany these visual displays, the portable 
computer also calculates and stores values for concentration path length with units of 
ppm-m which correspond to the real time displayed SO 2 absorption spectrum. 
Concurrently, the GPS collects position information and the computer stores these values 
in a file with the average concentration path length. This automated calculation of 
concentration path length significantly reduces the amount of time required and chance 
for inconsistencies in the data analysis segment. The information provided by GPS can 
be used to calculate the straight-line segments of the traverse without requiring a map. 
Figure 8 summarizes FLYSPEC’s Data Collection segment in terms of system inputs, 
components, and outputs. 
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(2) Data Analysis Segment. With the average concentration 
path length and straight-line segments determined in the data collection segment, wind 
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Input = \ 




m 

Solar radiation / 

UV Sensor 


Laptop Computer 




Global Positioning System 
with preprogrammed 
waypoints for “straight line” 
road segments 


Output = Real time UV Spectrum and 
S02 absorption profiles along with stored 
values for average S02 concentration 
path length and corresponding GPS 
coordinates for each reading 


Figure 8. FLYSPEC System Inputs, Components, and Outputs Used 

in the Data Collection Segment 


speed data is the only other variable needed to produce the average SO 2 emission rate 
with units of tons per day. The wind speed data is collected separately and then 
combined with the average concentration path length to produce the emission rate. 
Figure 9 illustrates the Data Analysis Segment in terms of system inputs, components, 
and outputs. The calculations required in the data collection and data analysis segments 
are discussed in further detail in Section d. 


Input = Average ' 
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Laptop Computer 
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Output = 
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Figure 9. FLYSPEC System Inputs, Components, and Outputs Used 

in the Data Analysis Segment 


b. Sensor Design 

FLYSPEC’s UV sensor is the USB2000 Miniature Fiber Optic 
Spectrometer developed by Ocean Optics, Inc, Dunedin, Florida. The company started in 
1989 when researchers in Florida developed a fiber optic pH sensor as part of an 
instrument designed to study the ocean’s role in global warming. As part of the sensor 
design, they wanted to place a spectrometer on a buoy. At the time, there were no 
spectrometers small enough to do this, so they developed their own. The end product 
was a fiber optic spectrometer nearly a thousand times smaller and ten times less 
expensive than existing systems. (Ocean Optics, p.l) The USB2000 is the second 
generation of these miniature spectrometers. Several physical and performance 
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specifications for the sensor are listed in Table 2. For comparison purposes, the 
following paragraphs will describe the sensor in the same three sections as were used 
with COSPEC: front optics, mid optics, and rear optics. Figure 10 displays the sensor 
components divided into these sections. 

(1) Front Optics. Solar radiation is collected through a filter 
lens that is built into the sensor’s carrying case. This filter aids in reducing the amount of 
stray light reaching the spectrometer. Normally, the radiation will then be collected with 
the sensor’s small telescope with a focal length of 42 mm. However, during calibration, 
the radiation will pass through the calibration cells before entering the telescope. 
FLYSPEC has a high and low calibration cell. The glass quartz cells have slightly 
different known concentrations than the cells used in COSPEC, but the same company 
manufactures them. To perform calibrations, the calibration cells are rotated into the path 
of the radiation in the same way demonstrated by COSPEC. The cells are mounted to the 




Figure 10. FLYSPEC Sensor Components (After Bo Galle, Figure 5, p.40) 


outside of the spectrometer and are rotated above the telescope as illustrated in Figure 11. 
Once the radiation has passed through the calibration cells and telescope, it passes 
through an entrance slit into the mid optics section. 

(2) Mid Optics. The radiation is directed and focused onto a 
diffraction grating that disperses the radiation and directs it to the rear optics. The 
components and functions of this section are very similar to COSPEC, but they operate 
on a much smaller scale. 
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(3) Rear Optics. The dispersed radiation is focused onto a 


High and Low 
Concentration 
Calibration Cells 


Figure 11. FLYSPEC Telescope and Calibration Cells 

2048-element charge coupled device (CCD). The CCD is a linear silicon array. It is 
covered in a coating that optimizes the CCD for detection of UV radiation. The CCD 
detects UV radiation and forms a spectrum that ranges from 177 - 333 nm. This entire 
spectrum is transferred to the portable computer. This is a significant difference from the 
data that COSPEC produces which corresponds to a limited number of specific 
wavelengths in the UV spectrum. A software program has been developed to display this 
detected UV spectrum, calculate and display a corresponding SO 2 absorption spectrum, 
and maintain a log of concentration path length values in ppm-m along with the 
corresponding GPS location recordings at a rate of 1 Hz. 

The following table presents a summary of physical and 


Physical Specifications 

FLYSPEC 

COSPEC 

Dimensions (LWH in inches) 

3.5 x 2.5 x 1.31 

31 x 14.74 x 8.75 

Weight 

.6 lb with cable 

42.5 lbs 

Power 

.5 W 

8 W @ 12VDC 

Operating Temp Range (Celcius) 

10 - 50 

0 - 50 


Detector Specifications 


Detector Type 

2048-element Linear Silcon CCD array 

Photomultiplier Tube 

CCD element size (microns) 

12.5 x 200 

N/A 

Effective Range (nm) 

177 -333 

280 - 320 


Optics Specifications 


Grating 

2400 lines 

1200 l/mm, plane, 68 x 54 mm 

Slit (WH in microns) 

25 x 1000 

.25x5.8 

Focal Length 

42 mm (input) 68 mm (output) 

25 cm 

Optical Resolution (nm) 

-.25 

0.4 

Pixel Resolution (nm) 

-.1 

N/A 

Stray light 

< 0.10% at 250 nm 

N/A 


Performance Specifications 


A/D resolution 

12-bit 

12-bit 

A/D sampling frequency 

1 Hz - 1 KHz 

1 Hz (max) 

Integration Time 

3 milliseconds to 60 seconds 

1 second - 32 seconds 

Field of View (degrees) 

10 

1.3178 x .401 

Sensitivity 

86 photons/count or 2.9*10 A -17 watts/count 

2.8 millivolts/ppm.m (typical low range) 

Signal to Noise Ratio 

250 : 1 -- single acquistion at full signal 

N/A 

S02 Detection Limit (low range) 

9-13 ppm.m 

< 5 ppm.m 

Table 2. Physical and Performance SpeciJ 

1 cations of the 


COSPEC and FLYSPEC sensors 
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performance specifications for the sensors on both COSPEC and FLYSPEC. There are 
several differences between the sensors, but some of the key differences are at the system 
level. The new spectrometer technology has offered solid improvements to SO 2 
detection, but it is the combination of new spectrometer technology and new computing 
technology that may allow FLYSPEC to prove the ability collect and process data more 
efficiently, consistently, and effectively than COSPEC. 

c. Data Collection Methods and Techniques 

FLYSPEC hopes to prove capable of collecting data using all three 
techniques demonstrated with COSPEC. To date, FLYSPEC has only demonstrated 
collection of data using the same ground mobile technique. FLYSPEC performs the 
ground mobile technique using a similar method to COSPEC; however, there are a few 
differences. Due to FLYSPEC’s smaller size, it is mounted to the top of the vehicle, 
either on the hood or on the roof, so that it can point directly up at the sky. FLYSPEC 
has made several traverses collecting data concurrently with COSPEC in this manner (see 
Figure 12). 

At the beginning of a traverse, while the sensor is under “clean air”, 
FLYSPEC collects several frames of data on each of its calibration cells just as COSPEC 
does. In addition to the calibration data, FLYSPEC also collects data on a dark frame and 
a reference frame. A dark frame is created by simply covering the lens of the carrying 
case so that no radiation can enter the sensor. A reference frame is created by looking 
straight up at “clean air”. All of this initial data collected prior to the start of the traverse 
plays a very important role in the data analysis and calculations. Following these initial 
measurements, a traverse is made. Calibration, dark, and reference data is occasionally 
collected at the end of a traverse; however this is not required. 



Figure 12. Concurrent Data Collection with COSPEC and FLYSPEC 
Using the Ground Mobile Technique 
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d. Data Analysis and Calculations 

Although the data collection methods and techniques of COSPEC and 
FLYSPEC are nearly identical, the differences in the technology of the sensors translate 
to slightly different data analysis and calculations. FLYSPEC calculates a real time SO 2 
absorption profile from data collected over a range of wavelengths using Equation 1.5. 
The calculation uses real time sample data, reference frame data, and dark frame data 
over the range of UV wavelengths to produce an absorption profile. 


A-x 


— — log 
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-AiJ 


(1.5) 


A x = Absorbance at some wavelength (unitless) 
S x = Sample intensity at some wavelength 


D x = Dark intensity at some wavelength 


R x = Reference intensity at some wavelength 
Figure 13 shows a UV absorption spectrum for a laboratory sample of SO 2 with several 



Figure 13. SO 2 Absorption Spectrum in UV Spectral Region 
(After Galle, Figure 2., p.17) 
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peaks and troughs annotated. Most of the strong SO 2 absorption peaks and troughs are 
located between 300 and 310 nm. Although FLYSPEC collects over the much broader 
range of 177-333 nm, many of the figures illustrating FLYSPEC data will be truncated 
to a more limited range of 300 - 310 nm in order to highlight these stronger absorption 
features of SO 2 . 

Figure 14 displays sample data from one of FLYSPEC’s traverses on 


FLYSPEC Reference Spectrum 





Figure 14. Example Reference, Dark, and Sample Data frames collected 
by FLYSPEC during traverse on Kilauea Volcano, Hawaii on March 4, 2002 
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Kilauea Volcano, Hawaii on March 4, 2002. Sample (a) illustrates a reference data frame 
during the traverse to produce real time absorbance spectra. Sample (c) illustrates a data 
frame in the middle of the traverse in which SO 2 is present. Although the absorbance 
spectra are calculated and displayed real time on the portable computer, these spectra are 


FLYPEC High Calibration Spectrum 




Chain of Craters Ed 3/4/02 Run 3 



Strands 

(c) FLYSPEC’s Concentration Path Length Vs Traverse Time 


Figure 15. FLYSPEC’s Calibration Cell Spectrum and Concentration 
Path Length Plot for a March 4, 2002, Traverse On Kilauea Volcano, Hawaii 
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not stored but simply used as data in the next step of calculating the concentration path 
length. These absorbance spectra are fitted using a non-linear quadratic to a laboratory 
SO 2 spectrum generated by the calibration cells. High and low calibration data collected 
during this March 4th traverse are displayed in Figure 15 as samples (a) and (b) 
respectively. This spectral fitting results in the average concentration path length of the 
observed plume in ppm-m. These calculations are performed at a rate of 1 Hz. The 
results are values of concentration path length over the time of the traverse. The 
concentration path length plot representing the March 4th traverse is also included in 
Figure 15. The point highlighted on the concentration path length plot corresponds to the 
absorbance calculated using the data frame in Figure 14, sample (c). In addition to the 
complete spectrum calculations, FLYSPEC can also be programmed to perform and store 
these same calculations at specific wavelengths. This is similar to what COSPEC does 
with the slits in its correlation disc; however, FLYSPEC’s peaks and troughs are not 
limited by a specific number of physical slits and can be easily changed through the 
computer’s custom built software. 

The concentration path length results are recorded along with 
corresponding GPS coordinates and are stored in a text format on the portable computer. 
All of these calculations are performed real time during the data collection. A sample of 
the typical FLYSPEC data text format is show in Figure 16. The time is Hawaii Standard 
Time. 


Data Point 

Year 

Month 

Day 

Hour 

Minute 

Second 

Latitude 

Longitude 

Elevation 

Raw S02 

Smooth S02 

0 

2002 

4 

18 

9 

58 

54 

1924.1201 N 

15515.9316 W 

1131.4 

5.93 

5.93 

1 

2002 

4 

18 

9 

58 

55 

1924.1201 N 

15515.9316 W 

1131.5 

5.8 

5.8 

2 

2002 

4 

18 

9 

58 

56 

1924.1202 N 

15515.9316 W 

1131.5 

5.2 

3.03 

3 

2002 

4 

18 

9 

58 

57 

1924.1202 N 

15515.9316 W 

1131.5 

-1.89 

-0.65 

4 

2002 

4 

18 

9 

58 

58 

1924.1202 N 

15515.9316 W 

1131.5 

-5.25 

-5.15 

5 

2002 

4 

18 

9 

58 

59 

1924.1202 N 

15515.9316 W 

1131.5 

-8.31 

-4.06 


Figure 16. FYSPEC Data Text Format Stored During Data Collection 


The latitude and logitude readings are in degrees and decimal minutes. The elevation 
readings are in meters. The raw and 3-point smooth readings of SO 2 concentration path 
length are in ppm-m. 

Following the data collection, FLYSPEC uses the same equation as 
COSPEC (Equation 1.3) to generate an emission rate in tons per day. The stored average 

concentration path length is combined with corresponding wind data and segment length 
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data. The wind data (speed and direction) is obtained from a ground station near the 
volcano’s vent. The segment length can be derived from the GPS data; however, the 
computer program is not fully developed yet and does not use the GPS data in this way. 
Currently a map is being used to generate the segment length which is the same method 
used in analyzing COSPEC’s data. 

These data collection and analysis methods are still under development. 
Several traverses have been made concurrently with COSPEC in order to validate the 
sensor and calculations. Figure 15 (c) is one of several comparisons that have been made 
between COSPEC and FLYSPEC data up to this date. In general, the results have been 
very similar and promising. There are a few minor differences that have not yet been 
fully explained. Some of differences are believed to stem from the substantial physical 
difference that exists between the fields of view of the two sensors. The variances in the 
results may actually correspond to the fact that each sensor is physically “seeing” a 
different “piece of sky”. Both sensors are also sensitive to different sun angles and cloud 
cover, and their degrees of sensitivity to these issues may be the cause of the different 
variances in the results. Additional tests and concurrent traverses with COSPEC will 
help eliminate many of these questions. 

B. SULFUR DIOXIDE EMISSION RATES OF KILAUEA VOLCANO, 

HAWAII, USING UV SENSORS 

Kilauea’s long-term activity has made it an ideal volcano for many different 
studies. Located in Hawaii, Kilauea currently hosts the longest running SO 2 emission- 
rate data set on the planet. (Sutton, Elias, Gerlach, and Stokes, 2001, p. 283) Stoiber and 
Malone first used COSPEC to measure Kilauea’s SO 2 emission rates in 1975. A program 
of nearly weekly SO 2 emission rates was started in 1979 and has continued through the 
present. (Casadvevall and others, 1987; Greenland and others, 1985; Elias and others, 
1993; Elias and Sutton, 1996). The US Geological Survey’s Hawaiian Volcano 
Observatory (HVO) has served as a headquarters for this SO 2 monitoring, and has 
published a compilation of COSPEC measured SO 2 emission rates covering the period 
from 1979 through 1997. Much of HVO’s staff and volunteers have participated in the 
monitoring program consisting of approximately 1100 days of measurements including 
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more than 5,000 plume traverses over the 19-year period. (Elias, Sutton, Stokes, 
Casadevall, 1998, p.3) 

Initially, in 1979, Kilauea’s main area of SO 2 release was its summit caldera. In 
1983, an eruption began along the East Rift Zone (ERZ) initiating an additional 
significant area of release at the Pu’u’O’o vent. Currently, these two areas are the main 
focus for SO 2 emission rate measurements. Measurements have been made in both of 
these areas using all three methods available to COSPEC; however, only the vehicle 
based traverse method results will be discussed, as this is the method of focus for this 
paper. 

Figure 17 illustrates the two main areas of SO 2 release on Kilauea. Measurements 
of SO 2 emission from the summit caldera are typically made along a section of Crater 
Rim Drive which encircles the caldera. This section of road is highlighted in Figure 17. 
Typically brisk northeasterly trade winds (3-12 m/s from 0 to 45 degrees east) create a 
plume of SO 2 that drifts approximately 1 km from its Halemaumau and solfatara source 
across Crater Rim Drive. Under ideal wind conditions, the plume will be dense, compact, 
and measurable close to the source. Daily emission rates are typically calculated from an 
average of data collected over six or more ten-minute plume traverses in the summit 
caldera area. 



Caldera and East Rift Zone (After Sutton et ah, Figure 1., p.284) 


23 













At this point it is important to note that wind speed is normally one of the largest sources 
of error in the data analysis. Several different methods have been used to measure wind 
speed over this extensive data collection period. Experiments have determined that wind 
measurements collected from 3 meters above the ground at a site next to HVO most 
accurately reflect the plume speed in the caldera region. As of 1998, wind speed 
measurements are made at a station 3 meters above ground at HYO. This station is 
labeled as the summit wind monitor site in Figure 17. Data collected in the past have 
been adjusted by up to 20% in some instances in order to normalize the variances that are 
believed to have been caused by the different wind collection methods. With the 
appropriate wind adjustments, Figure 18 displays SO 2 emission rates calculated for the 
Kilauea caldera area from June 1979 through 1997. The units of emission rates are 
metric tons per day (mt/d). An increase in emission rates begins with the eruption of 
Pu’u’O’o in 1983. A continuous decline is seen from 1987 through 1997. 
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Figure 18. Long Term Data Set Of Emission Rates Determined By 
COSPEC For The Kilauea Summit Caldera Area (From Elias et al., Figure 3., p.5) 


SO 2 measurements collected along sections of Chain of Craters road were 
initiated in 1992 in order to better account for the emission rate of Pu’u’O’o. These 
sections are also highlighted in Figure 17. Ideal wind conditions for these measurements 
consist of speeds greater than 5 m/s and directions between 25 and 40 degrees east. 
These conditions will normally create a compact plume approximately 9 km downwind 
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of the Pu’u’O’o vent which crosses the Chain of Craters Road above the 180-degree turn 
at the Holei Pali annotated in Figure 17. Daily emission rates for these sections are also 
calculated from an average of data collected over several plume traverses. 

HVO performed multiple experiments to determine the best method for collecting 
wind data in the Chain of Craters Road area. It was determined the winds measured near 
the center of the traverse and a good distance above the ground provided the best 
measurements. Currently, HVO uses a continuous wind monitor 3.5 m above the ground 
near the 180-degree turn at Holei Pali. This wind monitor site is also labeled in Figure 
17. Figure 18 displays emission rates calculated for the Chain of Craters Road section 
from 1992 through 1997. The black vertical bars represent the standard deviation of all 
traverses on a single day. The units for average daily emission rates are mt/d. (Sutton, 
Elias, Gerlach, and Stokes, 2001, p. 285-286) 
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Figure 19. Long Term Data Set Of Emission Rates Determined By 

COSPEC For The Chain of Craters Road Area (From Elias et al., Figure 4., p.6) 


C. SUMMARY OF UV DETECTION 


There are advantages and disadvantages in using COSPEC and FLYSPEC to 
monitor volcanic SO 2 emission rates. Both of the sensors operate using scattered 
skylight, which allows for simplistic and versatile collection methods and techniques for 
the sensors. Neither sensor requires a specialist operator or subsequent spectral analysis 
because the fit to the SO 2 spectrum is achieved with the calibration cells. The main 
disadvantage is that the emission rates are subject to large errors due to uncertainty in the 
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plume velocity. The speed of the plume is often assumed to be equal to wind-speed 
measurements obtained close to the ground. These wind speeds near the ground are not 
always accurate for the plume. (McGonigle and Oppenheimer, 2002, pp. 4-5) 

The concepts used for the data collection and analysis are very similar for 
COSPEC and FLYSPEC; however, the advantages offered by the new technologies used 
in FLYSPEC should not be overlooked. The most obvious advantage is in the physical 
size of the sensor. COSPEC is roughly the size and weight of a 40-pound sack of 
potatoes, whereas FLYSPEC is roughly the size and weight of a deck of cards. Some 
volcanic plumes are accessible only on foot. Having a sensor the size of FLYSPEC to 
complete this form of data collection and analysis would greatly simplify the process. 
FLYSPEC is also less expensive which could facilitate the ability of researchers to obtain 
multiple sensors. COSPEC collects data corresponding to a limited number of specific 
wavelengths within the UV spectrum, whereas FLYSPEC collects data over an entire 
band of the UV spectrum. This range of data combined with the spectral fitting algorithm 
could prove to provide a more accurate and consistent means of analyzing data. The 
addition of the new computing power to the system makes this analysis possible in real¬ 
time. The lower cost, ability to have multiple sensors, automatic analysis, and data 
storage capability could facilitate putting several FLYSPECs at a remote site to collect 
SO 2 plume samples over several days at the same spot. These are significant advantages 
which will hopefully be proven over time and provide researchers a better, faster, 
cheaper, and smaller means of detecting and analyzing volcanic SO 2 emission rates. 
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III. INFRARED DETECTION OF S0 2 

A. AIRBORNE AND SPACEBORNE INFRARED SENSORS 

Although COSPEC has set the precedence for detecting and measuring 
concentrations of S0 2 within volcanic plumes, its various data collection and analysis 
techniques do have some inherent limitations. Most of these limitations originate in the 
geographic locations of the various volcanoes, in the volcanoes’ natural and man-made 
surroundings, or in the standard meteorological conditions of the area, rather than in the 
COSPEC system or sensor. For example, if typical wind conditions do not generate a 
plume drift with appropriate geometries to an existing road system, COSPEC’s ground 
mobile collection technique could become difficult and produce less accurate results. 
COSPEC’s airborne collection technique requires the sensor to be flown beneath the 
plume, and if the plume consistently drifts near the ground, an aircraft may not be able to 
achieve an altitude low enough to allow COSPEC to collect data. Solutions to these 
limitations could be developed through new data collection and analysis techniques; 
however, they may also be discovered in the application of sensors operating in a 
different region of the electromagnetic spectrum, such as thermal infrared sensors. 

Collecting data in a different region of the electromagnetic spectrum not only 
implies different data collection and analysis techniques, but it also generates the 
possibility of seeing new characteristics in the data. In addition to measuring the plume’s 
S0 2 concentration, the transition to the thermal infrared region gives the sensor the ability 
to collect information on thermal aspects of the volcano. Several airborne and 
spacebome infrared sensors have recently been used to explore some of these 
capabilities. 

NASA has used its airborne Thermal Infrared Multispectral Scanner (TIMS) to 
develop data analysis techniques and potential algorithms to map volcanic plumes and to 
determine a plume’s S0 2 concentration and emission rate. TIMS has also been used to 
map temperature and emissivity of volcanic lava surface flow fields. These studies with 
TIMS laid initial groundwork that is being further developed for use with the Advanced 
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Spacebome Thermal Emission and Reflection Radiometer (ASTER), a sensor currently 
flown on a Terra Satellite that is part of NASA’s Earth Observing System (EOS). 

The Airborne Hyperspectral Imager (AHI), owned by the University of Hawaii, is 
another thermal infrared (TIR) sensor that can be used to monitor volcanic activity. This 
paper explores AHI’s ability to detect volcanic plumes and determine a plume’s SO 2 
concentration. The following paragraphs provide an overview of the TIMS sensor, its 
data analysis techniques, and some of its published results. ASTER and its potential to 
use the TIMS analysis techniques are also discussed. Finally, a description of AHI’s 
system and sensor design, along with its data collection techniques is presented. 

1. Thermal Infrared Multispectral Scanner (TIMS) 

TIMS is an airborne long wave infrared (LWIR) sensor developed by NASA 
Stennis Space Center, the Jet Propulsion Laboratory (JPL), and Daedalus Corporation. 
TIMS has demonstrated abilities in a variety of applications from ground temperature 
mapping and mineral classification to detecting volcanic plumes and determining SO 2 
concentrations. In addition to performing its airborne remote sensing mission, TIMS has 
also served as a simulator for ASTER, a spacebome IR sensor. A brief overview of the 
sensor follows, as well as a description of data analysis techniques that allow the sensor 
to detect and quantify volcanic SO 2 plumes. 

a. Sensor Design and Data Collection 

TIMS is a multispectral LWIR sensor with a six-element HgCdTe array. 
This array is sensitive to an LWIR band ranging from 8.2 to 12.2 microns. A dispersive 
grating spreads this LWIR band over the six pixels so that each pixel is sensitive to a 
smaller segment of the broad LWIR region. (Baer-Riedhart, “TIMS”) 

TIMS operates from several different aircraft: C-130, ER-2, Stennis 
Learjet, and DOE Cessna Citation aircraft. It operates in a scanning mode and scans the 
array over a 76.56 degree FOV collecting data for 638 points on the ground per scan line. 
It has a variable scan rate from 7.3 to 25 scans per second. (Hydrology Data Support 
Team, “Performance Parameters”) The multispectral image cube that is produced 
consists of 638 pixels along the x-axis, a number of pixels in the y-axis that are 
determined by the number of scan lines collected along the flight line, and six x-y planes 
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along the z-axis that each represent one of the smaller segments of the 8.2 - 12.2 
wavelength region. Figure 20 illustrates the scanning collection mode. Often in a 



6-demat FPA 



Iniqy length (ktemined Iw 


Ground 


the rubber of scan lints collected 



# of pixels in y-axis determined 
by number of scan lines 



Figure 20. Scanning Multispectral Data Collection And 

Resulting Data Cube 


scanning mode the array is stationary, and some form of rotating mirror performs the 
actual scanning motion; however, the figure depicts TIMS’ six-element array performing 
the scanning motion. This was done in order to clearly illustrate the relationship between 
the scanning array and resulting multispectral data cube. 


Detector Specifications 


Imager 

6 element HgCdTe FPA 

Spectral Coverage (microns) 

8.2-12.2 

Number of Spectral Bands 

6 


Optics Specifications 


Spectral Resolution 

700 nm 

Angular Resolution (mrad) 

2.5 

Ground Resolution (m) 

50 at 65,000 ft altitude 

Aperture (in) 

7.5 


Performance Specifications 


Quantization 

8 bits 

Pixels per scan line 

638 

Scan Rate (scans/second) 

varable: 7.3, 8.7, 12, 25 

Software 

Data can be read with any standard 

image processing package, e.g. PCI, ERDAS, IDL etc. 

File Size 

4-9 MB 

Set Size 

-2.5 GB 

FOV (full / angular) 

76.5 degrees 

Swath Width (km) 

31.3 

Guidance Sensor 

roll corrections of + / -15 degrees 

Operating Alitudes (km) 

1 -20 

NEDT 

< .4 K 

Table 3. TIM 

IS Sensor Specifications 
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The sensor has a Noise Equivalent Delta Temperature (NEDT) of less than 
0.4 degrees K. Two on-board black bodies perform calibration for the sensor. High and 
low temperatures are chosen for these black bodies based upon the temperatures expected 
to be imaged in a particular scene. Data is collected on these black bodies at the 
beginning and end of each scan line. Additional sensor specifications are included in 
Table 3. (Baer-Riedhart, “Sensor/Aircraft Parameters”) 
b. Data Analysis Techniques 

As seen with FLYSPEC, data analysis of SO 2 plume concentrations using 
ultraviolet wavelengths is based on an absorbance profile calculation over a straight-line 
distance transecting the plume. These results are fitted to a known laboratory SO 2 
concentration absorbance profile. The data is collected looking up at the sky, providing a 
fairly constant background to measure against, with the main exception being cloud 
cover. Data analysis of SO 2 plume concentrations using thermal infrared wavelengths is 
slightly more complicated than using ultraviolet wavelengths. Calculations are based on 
a form of the radiative transfer equation. Using TIMS or AHI, the data is collected 
looking down at the ground, possibly resulting in measurements against a highly varied 
background. 

(1) Estimating SO 2 Concentration. Realmuto et al. (1994) 
introduced a procedure for estimating the SO 2 content of volcanic plumes using data 
collected by TIMS. This procedure models the radiance observed by TIMS as it looks to 
the ground through an SO 2 plume using the MODTRAN radiative transfer code (Berk et 
al., 1989). 

MODTRAN calculates atmospheric transmission and radiance using 
models of the absorption bands of 12 gas molecules (H 2 O, CO 2 , O 3 , N 2 O, 

CO, CH 4 , O 2 , NO, SO 2 , NO 2 , NH 3 and HNO 3 ), utilizing three-parameter 
band models with spectral resolutions of 1 cm' 1 . Realmuto et al., 1997, 
p.15059) 

Equation 3.1 is the basic radiance equation used in the code. It is a simplified version of 
the equation, as it does not account for any viewing angles other than nadir. 

L S {A,T 0 ) = {e{A)B{h,T,)+[\ - e{X)]L d {X)W) + L U {X) ( 3 . 1 ) 
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L s {A, T 0 ) = Radiance of the ground as it is viewed through the atmosphere as a 

function of wavelength (A ) and ground temperature (To) 

e(A) = Ground emissivity as a function of wavelength 

B(A, T 0 ) = Planck function in terms of wavelength and ground temperature 

L d {X) = Ambient or sky radiance at ground altitude as a function of wavelength 

t(A) = Spectral transmittance of the atmosphere as a function of wavelength 

L u ( A ) = Ambient or sky radiance at instrument or sensor altitude as a function of 
wavelength 

Table 4 lists inputs and their sources that are required to run the MODTRAN code. With 
these inputs, the model calculates values for L s , Ld, Lu, and T . Figure 11 illustrates two 
example atmospheric transmission profiles generated by MODTRAN. The dashed line 
represents a transmission profile through an atmosphere that is free of SO 2 whereas the 
solid line represents a transmission profile through an atmosphere that contains SO 2 . A 
typical SO 2 plume is cooler than the ground so it will absorb the ground radiance making 


Input Parameter 

Source 

Ground emissivity 

Estimate from TIMS radiance measurements of ground that is not beneath the plume 

Ground altitude 

DEM or other mapping product 

TIMS or sensor Altitude 

Aircraft measurements 

Zenith angle 

Path between sensor and ground determined by aircraft measurements 

Barometric pressure 

Meteorlogical measurement or packaged LOWTRAN atmospheric model 

Temperature 

Meteorlogical measurement or packaged LOWTRAN atmospheric model 

Humidity 

Meteorlogical measurement or packaged LOWTRAN atmospheric model 


Table 4. MODTRAN’s Required Inputs And Respective Sources 


it an identifiable feature in the TIMS ground radiance profiles. This is illustrated by 
comparing the dashed line profile (0 gm" 2 of SO 2 ) and the solid line profile (23 gnf 2 of 
SO 2 ) in Figure 21. 


Indicated by shaded regions, Figure 21 also illustrates the 
normalized spectral responses for the six TIMS channels operating in a typical 
atmosphere free of SO 2 . The figure illustrates that if there is SO 2 in the atmosphere, the 
strongest response is featured in TIMS’ channel 2. Having the strong absorption 
response in one channel generates opportunities to expose the feature through applying 
different processing techniques that enhance or emphasize that particular channel. 
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Figure 21. Atmospheric Transmission Profiles With and Without SO 2 
(from Realmuto et al., Figure 2, p.483) 

Often multispectral datasets exhibit high correlation and produce 
rather bland color images; however, a decorrelation stretch can be applied to produce a 
more colorful composite image, which could more effectively display information 
contained in the image by applying colors to specific bands of wavelengths. A 
decorrelation stretch requires three bands for input. The bands will be represented in the 
red, green, and blue (RGB) colors. (Research Systems, 2001, p. 580) With TIMS data, 
bands 5,3, and 2 are often chosen to represent red, green, and blue respectively. The SO 2 
will absorb the ground radiance in channel 2, implying that the blue color is absorbed in 
the image. This will create an image with the SO 2 plume displayed in shades of yellow 
and red. This decorrelation stretch technique is illustrated in Figure 22 with a TIMS 
image acquired over the East Rift Zone of Kilauea Volcano, Hawaii on September 30, 
1988. (Realmuto et al., 1994, p. 482) 

The SO 2 absorption feature is a function of all of the input 
parameters listed previously in Table 4. These parameters can all be determined through 
various ancillary sources. There are two additional parameters that influence the 
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Figure 22. Decorrelation Stretch Applied To TIMS Image Data 
(Realmuto et al., Figure 2, p. 15060) 

existence of the SO 2 absorption feature that cannot be determined by ancillary sources: 
ground temperature and SO 2 concentration. SO 2 concentration is modeled as an 
atmospheric constituent in the MODTRAN code, and the ground temperature is estimated 
using the collected radiance data. There are many combinations of ground temperature 
and SO 2 concentration that will produce the same fit to a radiance spectrum, so the 
problem of quantifying the SO 2 absorption feature remains underdetermined at this point. 
Realmuto et al. (1994) uses a two-step process to produce a solution. (Realmuto et al., 
1997, p. 15061) 

The first step involves using TIMS radiance data from channels 4, 
5, and 6 to determine values for ground temperature. As is shown in Figure 21, these 
three channels are not affected by the presence of SO 2 in the atmosphere. Their 
“resilience” to the presence of SO 2 implies that the ground temperatures determined 
within these wavelengths should not fluctuate whether a plume is present or not. Several 
temperatures are estimated and used in the MODTRAN code with wavelengths 
corresponding to channels 4, 5, and 6. The SO 2 concentration is modeled as 0 gm~ 3 . The 
various MODTRAN radiance profile results are compared with the TIMS radiance 
profiles observed in channels 4, 5, and 6 until the model ground temperature that 
produces the best weighted least squares fit between the two radiance profiles is 
established. This model ground temperature is then used in the second step of the 
solution. (Realmuto et al., 1997, p. 15061) 
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The second step involves using TIMS radiance data from channels 
1, 2, and 3 to determine SO 2 concentration. The model ground temperature from the first 
step is used along with several estimated SO 2 concentrations. These are used in the 
MODTRAN code, and the radiance profile results are again compared with the TIMS 
radiance profiles. The results are compared until the SO 2 concentration that produces the 
best weighted least squares fit between the two is established. Channels 1 and 6, located 
near water vapor and CO 2 absorption bands respectively, are given less weight in both of 
the fitting processes. The SO 2 concentration determined in the second step can be 
multiplied by the plume thickness to produce an estimated SO 2 column abundance with 
units of gm' 2 . (Realmuto et al., 1997, p.15061) 

In addition to applying this estimation procedure to single pixels, it 
can also be applied to a line of pixels that transects the plume. This is similar to the way 
COSPEC collects data in straight-line distances perpendicular to the plume. Following 
the estimate of the column abundances for each pixel, the total SO 2 burden (gm' 1 ) can be 
integrated across the line of pixels for each transect. Multiplying the SO 2 burden by the 
wind velocity produces an estimated emission rate. An average of the estimated emission 
rates for each transect produces an emission rate for the day similar to the methods of 
COSPEC. (Realmuto et al., 1994, p. 485) 

(2) Plume Mapping Procedure. In 1997, Realmuto et al. extended 
the SO 2 estimation procedure to develop a plume mapping procedure. The estimation 
procedure is based on data from individual pixels and the mapping procedure simply 
extends this to allow an estimation of SO 2 over blocks of adjacent pixels, creating a two- 
dimensional SO 2 plume map. 

Realmuto et al. applied this mapping procedure to data acquired in 
an area of the East Rift Zone of Kilauea Volcano on September 30, 1988. Pu‘u‘0‘o was 
one of the active vents providing data in 1988, and it continues to be active today. Figure 
23 illustrates a map of the ground temperature estimation step with a corresponding SO 2 
column abundance estimation map for the Pu’u’O’o vent. The quantified results of SO 2 
burden and emission rate for each transect of the plume are displayed in Table 5. 

Realmuto et al., draws several conclusions about the scene by 
comparing the ground temperature and column abundance maps. The ground 
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temperature map indicates variations in the ground temperature beneath what appears to 
be a fairly constant SO 2 column abundance on the west edge of the plume. This anomaly 



Figure 23. Estimation Maps of Ground Temperature and SO 2 Column 
Abundance Generated For Pu’u’O’o Vent Of Kilauea Volcano, Hawaii 
(after Realmuto et al., Figure 3 and Figure 4, p. 15061 - 15062) 


T ransect 

Total S02 Burden 

g/m 

S02 Emission Rate 

kg/s t/d 

A 

5176 

29.5 

2549 

B 

3763 

21.4 

1849 

C 

1989 

11.3 

976 

D 

2505 

14.3 

1236 


Table 5. Quantified TIMS Results Of Pu’u’O’o Data Collected 
On September 30, 1988 (After Realmuto et al., 1997, p. 15066) 


was determined to be caused by the sensitivity of the ground temperature estimates to the 
different textures of lava in the scene. The rough surface of an aa flow allows 
aerodynamic cooling so that it is cooler than glassy pahoehoe flows in the ground 
temperature maps. Figure 24 displays these two common types of volcanic surface 
textures. (Realmuto et al., 1997, p. 15064) 



Figure 24. Typical Pahoehoe And Aa Lava Flows 


The SO 2 abundance map indicates that there are small areas of the 
plume that were significantly denser than others. These areas are referred to as puffs 
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within the plume. A comparison between these puffs and the corresponding ground 
temperatures produces several anomalies indicated on the ground temperature map. 

The northern puffs do not create the expected variation in ground 
temperature, although a small variation is present. It was determined that these northern 
ground temperature anomalies are more likely to have been produced by a small amount 
of water vapor in the puffs that was not accounted for in the estimation procedure. These 
ground temperature variations are not necessarily related to the apparent increase of SO 2 
column abundance in the puffs. There is not sufficient evidence to invalidate either the 
ground temperature or column abundance results, so they are maintained and used in the 
final emission rate calculations for the vent. (Realmuto et al., 1997, p.15064-15065) 

The southern puffs do appear to cause a significant variation in 
ground temperature. With the aid of data collected in visible wavelengths, it was 
determined that these puffs actually correspond to meteorological clouds that were 
present at the time of data collection. These puffs are therefore excluded from the final 
emission rate estimate. These examples stress the importance of comparing temperature 
and colu mn abundance results in order to detect and validate any anomalies. (Realmuto et 
al., 1997, p. 15065) 

(3) Comparison of TIMS Plume Mapping Procedure With 
COSPEC Results. The HVO Staff collected COSPEC data for the Pu’u’O’o vent using a 
tripod and the stationary ground collection technique on September 20, 1988, just 10 days 
before the collection of TIMS data used in producing Figure 23 and Table 5. COSPEC 
collected data for 29 transects over 3.5 hours, while TIMS collected data for 4 transects 
over 16 seconds. This difference in time required for data collections can have 
significant impacts on the measurement results as is discussed in a following paragraph. 
Table 6 is a comparison of these TIMS and COSPEC results collected over the Pu’u’O’o 
vent. (Realmuto et al., 1997, p. 15067-15069) 

Although the emission rate is the standard reported COSPEC 
measurement, Realmuto et al. (1997) state that a comparison between the results of TIMS 
and COSPEC would be more appropriate using the SO 2 burdens of individual transects. 
Several reasons are given to support this statement. For both instruments, the emission 
rate is not directly measured. It is the result of multiplying the SO 2 burdens of the 
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Measurement Units 

kg/s 

t/d 

Average Daily Emission Rate 



TIMS (Sept 30, 1988) 

16.8-20.2 

1449- 1741 

COSPEC (Sept 20, 1988) 

7.6-11.6 

661 - 1001 

Average Yearly Emission Rate 



COSPEC (1988) 

8.3-12.9 

720- 1120 

Highest Emssion Rate For Year 



COSPC (1988) 

14.5 

1255 


Table 6. Comparison of TIMS and COSPEC Results 


transects by wind speed. As noted in an earlier chapter, wind speed is a large source of 
error in the COSPEC measurements. The SO 2 burdens of COSPEC’s transects ranged 
from 1100 to 3000 gm-1. The TIMS SO 2 burdens ranged from 1900 to 5100 as seen in 
Table 5. This is a much closer comparison than that shown by the emission rates in Table 
6, promoting the comparisons of SO 2 burdens vice emission rates. 

Realmuto et al. (1997) suggest that the small difference that exists 
in the SO 2 burdens could stem from the puffs within the plume and the collection 
techniques of both instruments. The slow speed of COSPEC’s ground stationary 
collection method and the low wind speeds might have allowed the puffs to be diluted 
within the measured profile. For TIMS collection, wind speed in the direction of flight 
was nearly negligible relative to the ground speed of the aircraft. This meant that the 
puffs would have appeared instantaneous to TIMS accounting for the higher SO 2 burdens 
in TIMS’ transects A and B. 

(4) Errors and Sensitivity of Estimation and Mapping Procedures. 
In 1997, Realmuto et al. determined an error budget of +/- 20% for the mapping 
procedure. This is considered a minimum error, as not all of the error sources can be 
quantified. The main sources of error noted by Realmuto et al. are: the accuracy of 
radiance measurements made by TIMS, approximations such as plume thickness and 
altitude that are made during the estimation procedure, and the description of the transfer 
of radiation from the ground to the sensor involving atmospheric conditions and ground 
emissivity. (Realmuto et al., 1997, p. 15067) 

Several analyses have been conducted to determine the sensitivity 

of the mapping procedure to the various error sources. It is noted that the procedure is 

most sensitive to estimates in the plume altitude. Realmuto et al., (2000) suggest several 

methods to determine plume altitude when direct field observations cannot be made. An 
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assumption that the plume altitude is approximately equal to the elevation of the source 
vent is one of the simplest methods. Shadows cast by the plume may be used with the 
solar azimuth and elevation at the time of data collection as another method. Glaze et al. 
(1999) have developed a method based on photoclinometry, using the brightness of pixels 
to determine orientation of the corresponding surface to the sun. Regardless of the 
method used, the mapping procedure assumes that the plume is in thermal equilibrium 
with the troposphere and therefore defines the temperature contrast between the SO 2 
plume and ground. As the altitude estimates increase, there is a corresponding increase in 
temperature contrast. The increased temperature contrast decreases the amount of SO 2 
required to indicate absorption. Analysis determined that variances in plume altitude 
could cause errors as large as 12%. Variances in other sources caused errors on the order 
of 4% or less. Table 7 lists several error sources with corresponding results generated 
from various sensitivity analyses. (Realmuto et al., 1997, pp. 15066-15067) 


Error Source 

Variances Tested 

Sensitivity Results 

Plume Dimensions 



Plume Altitude 

Variances of up to 50% 

Decreases in S02 column abundnaces up to 12% 

Plume Thickness 

Reductions of up to 60% 

Decreases in S02 column abundnace less than 4% 

Atmospheric Profiles 



Relative Humidity 

Increase of up to 100% 

Decreases in S02 column abundance less than 4.3% 

Relative Humidity 

Reduction of 50% 

No change in S02 column abundance 

Model Atmospheres 

Substituted MODTRAN model atmosphere for radiosonde data 

Changes in S02 column abundance less than 10% 

Ground Emissivity 

Assumed constant emissivity 

Increases in S02 column abundances less than 4% 

Ground Elevation 

Reductions of up to 24% 

Decreases in S02 column abundance less than 1.5% 


Table 7. Sensitivity Analysis Results For SO 2 For Plume Mapping Procedure 


2. Advanced Spaceborne Thermal Emission Reflection Radiometer 
(ASTER) 

In 2000, Realmuto et al., extended the application of the SO 2 mapping procedure 
to data simulating satellite collection. For the collection simulation parameters, 
Realmuto chose to model the Moderate Resolution Imaging Spectroradiometer (MODIS) 
and Advanced Spaceborne Thermal Emission Reflection Radiometer (ASTER). These 
are two LWIR sensors employed on the Terra satellite that is part of NASA’s Earth 
Observing System. ASTER has a spatial resolution of 90 meters, much finer than the 1 
km resolution of MODIS. This paper will focus on the ASTER sensor, as its enhanced 
spatial resolution provides better comparisons between airborne and spaceborne sensors. 
a. ASTER Description 

ASTER is a cooperative effort between NASA and Japan’s Ministry of 
Economy Trade and Industry (METI). It is an imaging sensor with a main mission of 
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collecting data for detailed maps of land surface temperature, emissivity, reflectance, and 
elevation. ASTER has a Visible and Near Infrared (VNIR) subsystem, a Shortwave 
Infrared (SWIR) subsystem, and a Thermal Infrared (TIR) subsystem. Each subsystem 
has its own telescope. Table 8 displays characteristics of each of these subsystems. The 
focus for this paper is the TIR subsystem. 


Characteristic 

VNIR 

SWIR 

TIR 

Spectral Range (microns) 

.52 - .86 

1.6-2.36 

8.125-11.65 

Spectral Bands 

3 

6 

5 

Ground Resolution (m) 

15 

30 

90 

Data Rate (Mbps) 

62 

23 

4.2 

Cross-track Pointing (degrees) 

+/- 24 

+/- 8.55 

+/- 8.55 

Cross-track Pointing (km) 

+/-318 

+/-116 

+/-116 

Swath Width (km) 

60 

60 

60 

Detector Type 

Si 

PtSi-Si 

HgCdTe 

Quantization (bits) 

8 

8 

12 


Table 8. Characteristics of ASTER’s Three Subsystems 


(after Hook, “Instrument”) 


The TIR subsystem is sensitive to a wavelength range of 8 - 11 microns. 
This range is spread over 5 spectral bands. The system response of each of these bands is 
illustrated in Figure 25. Each band has 10 HgCdTe detectors in a staggered array. These 


Protoflight Model - System Response Functions 



- Channel 10 - Channel 11 - Channel 12 

- Channel 13 - Channel 14 


Figure 25. Spectral Response for ASTER’s TIR Subsystem 
(from Hook, “Instrument”) 

detectors are maintained at a temperature of 80 degrees K by a split Sterling mechanical 
cooler. The NEAT is less than .3 K for all bands. The sensor is calibrated with an 
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internal adjustable temperature black body. The TIR subsystem is fixed and requires a 
mirror to perform any scanning or pointing. 

b. ASTER Simulation and Results 

In the spaceborne simulation of the estimation and mapping procedures, 
two main differences emerge by moving the sensor’s operating altitudes from 1-20 km 
(TIMS) to 705 km (ASTER): the sensor’s resolution is decreased and the longer path of 
travel through the atmosphere decreases the amount of radiance perceived by the sensor. 
Realmuto et al. chose to work with the same TIMS, September 30, 1988, Pu’u’O’o data 
discussed previously. In order to create an image as it might be viewed through ASTER, 
the TIMS data was resampled to a 90-meter spatial resolution. After the data was 
resampled, the estimation and mapping procedures were applied. Figure 26 is an image 
of the area used in the simulation with a decorrelation stretch applied. As previously 
discussed, this causes the SO 2 plume to appear in shades of yellow and red. 



Figure 26. Simulation Data With A Decorrelation Stretch Applied 
(after Realmuto et al., Plate 1, p. 106) 


The SO 2 abundance results of the simulation are displayed in Figure 27. ASTER proved 
to show favorable results when compared with TIMS. The decrease in ASTER’s spatial 
resolution reduces or dilutes the column abundance estimates by increasing the area 
covered by each pixel. Although this effect did prohibit ASTER from producing some of 
the higher column abundances of the puffs that can be seen with TIMS, the dilution was 
not significant enough to cause the presence of the puffs within the plume to be 
completely overlooked. (Realmuto et al., 2000, p. 108) 
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Figure 27. Comparison of Pu’u’O’o Plume Maps At Spatial Resolutions 
of TIMS and ASTER (After Realmuto et al., Plate 1, p.106) 


With the concern of spatial resolution set aside, Realmuto et al. next 
addressed the concern of ASTER’s sensitivity. Another simulation was run to determine 
the apparent reductions in ground temperature that a sensor operating at 705 km might 
detect by viewing the ground through the 1988 Pu’u’O’o plume. The simulated 
transmission path was generated through a combination of radiosonde data (0-10 km) 
and the MODTRAN tropical atmosphere model (10 - 705 km). The simulation results 
are displayed in Figure 28. The changes in temperature are derived from radiance spectra 


Figure 28. 


Pu'u 'O'o Sensitivity Simulation 



Results of ASTER’s Sensitivity Simulation (After Realmuto et al., 
Plate 1, p.106) 


calculated for an atmosphere that was free of SO 2 and for three additional atmospheres 
with different simulated plume abundances. As previously mentioned, ASTER’s NEAT 
is less than .3 Kelvin for each of its 5 spectral bands. The temperature differences 
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derived for the 10 gm' and 15 gm' plume abundances exceed this value and would be 
detected; however, the smaller 5 gm' plume abundance does not generate a large enough 
temperature difference for ASTER to detect. (Realmuto et al., 2000, p. 108) 

With the promising results of the simulation, Realmuto et al. concluded 
that ASTER should be able to detect SO 2 plumes of similar size and characteristics to 
Pu’u’O’o. Realmuto et al. (2000) also suggest that ASTER could be used in conjunction 
with MODIS and TOMS (an ultraviolet spacebome sensor) for many important large 
scale studies such as: establishing baseline emission rates and performing change 
detection for volcanoes around the world, mapping atmospheric products generated by 
volcanic eruptions, and documenting the life cycle of SO 2 through its development to 
H 2 SO 4 . (Realmuto et al., 2000, p. 112) 

3. Airborne Hyperspectral Imager (AHI) 

The University of Hawaii’s Hawaii Institute of Geophysics and Planetology 
(HIGP) designed and built AHI as part of the Defense Advanced Research Projects 
Agency’s (DARPA) Hyperspectral Mine Detection (HMD) program. The main goal of 
the program was to develop and demonstrate a hyperspectral infrared capability for 
remote buried mine detection. The program outlined several high level system and 
sensor requirements based on an extensive phenomenology measurement program that 
guided the design and construction of AHI. (Lucey et al., 1998, p. 36) 
a. System Design 

The HMD program established the following five high-level system 
design requirements (Lucey et al., 1998, p. 38) 

• Must be airborne 

• Must provide real-time detection of mines 

• Must be demonstrated in a realistic scenario 

• Must demonstrate a clear path to an operational system 

• Must have the capability of serving as a phenomenology data collection 
platform 
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These requirements led to the development of the AHI system pictured in Figure 29. The 
system collects data through an infrared and visible sensor that are coupled with a real 
time processing system developed by Technical Research Associates, Inc. (TRA) and 



Figure 29. AHI System Components (after Lucey et al., Figure 3, p. 41) 


Space Computer Corp. The processing system controls many important system 
functions. It performs sensor control, generation of calibration coefficients, geometric 
preprocessing, real time radiometric calibration, near real time detection algorithm 
processing, user interface control, and data recording. A Xilinx field programmable gate 
array (FPGA) and four Share digital signal processors (DSP) are the main components of 
the processing system. These components receive the raw 12-bit digital data from the 
sensors, spectrally bin the data, apply calibration components, compute principal 
components, apply the detection algorithm, and then pass the calibrated raw data and 
principal components to a Pentium PC and 12 Gigabyte Redundant Array of Independent 
Disks (RAID) for storage. A GPS receiver also collects and sends information to the 
RAID for storage. The visible and infrared data can also be displayed as it is collected in 




















































waterfall plots on the system’s monitor. Figure 30 displays the system as it is mounted in 
a twin otter fixed wing aircraft while Figure 31 is a snapshot of the monitor display 
during data collection. (Williams “Sensor Details”) 



Figure 31. Snapshot of AHI’s Monitor Display 


b. Sensor Design 

In addition to system requirements, DARPA’s HMD program also 
outlined the following sensor requirements (Lucey et al., 1998, p. 38): 

• Spectral Range: 8.3 - 11 microns 

• Spectral Resolution: <32 wavenumbers 

• Noise Equivalent Spectral Radiance (NESR): < 0.01 watt/m 2 -micron-sr at 
300 degrees Kelvin 

• Pixels per target: >30 

With these baseline requirements, the University of Hawaii’s HIGP built a 
sensor that consists of the following four subsystems illustrated in Figure 32: telescope, 
spectrograph, background suppressor, and FPA with associated electronics. The 
telescope is a two-element diffraction limited transmission lens with a 111-micron focal 
length and 35 mm aperture. The spectrograph is an uncooled commercial reflective f/4 
imaging spectrograph with gold-coated optics. It is a grating spectrograph that has a 
linear dispersion at the output. An uncooled spectrograph is possible due to the sensor’s 
background suppressor. The background suppressor cools the FPA, a three-element 
transmission reimaging system, and a linear variable filter that are all housed in a vacuum 
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dewar. The FPA is a 256 by 256 Rockwell TCM2250 HgCdTe array sensitive to an 
LWIR band from 7.5 to 11.5 microns. It is cooled to 56 degrees Kelvin (K) by a 1.5 watt 



Figure 32. Main Subsystems of the Sensor Design (from Lucey et al., 

Figure 1, p. 37) 

split Stirling mechanical cryocooler. At 56 K, the half power drop-off in sensitivity is at 
11.5 microns. The array is operated at its maximum frame rate of 150 Hz with an 
integration time of 3 milliseconds in order to avoid saturation for scenes that may contain 
temperatures of 70 degrees Celsius. A separate 0.1 liter liquid nitrogen reservoir cools 
the sensor’s optics. (Lucey et al., 1998, p. 40) 
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Figure 33. Main Components of Sensor Subsystems (after Lucey et al., 
Figure 3, p. 41) 


The sensor also has an on-board adjustable temperature black body that 
provides information for real-time data calibration and conversion to radiance. In 


45 








































addition to the LWIR capabilities, the sensor can also collect data in visible wavelengths 
through a 3-color CCD linescan camera. These major components of the sensor’s 
subsystems are illustrated in the block diagram in Figure 33. (Lucey et al., 1998, p. 42) 

AHFs on-board computer has custom software that controls and monitors 
several sensor functions such as: control of the pod environmental shutter, the black 


Physical Specifications 


Size (LWH in inches of imager only) 

47 x 27 x 17 

Imager Weight (pounds) 

300 

Control & Rec Module (pounds) 

122 

Power Requirements 

1000 W (35 Amps @ 28 VDC) 

Operating Environment 

0 - 15,000 ft pressure altitude, -20 to +50 degrees C, 

Passive vibration isolation 

FPA Cooler 

Closed cycle sterling cooler - 56K 


Detector Specifications 


Imager 

256 x 256 element Rockwell TCM2250 HgCdTe FPA 

Spectral Coverage (microns) 

7.5-11.5 

Number of Spectral Bands 

256 or 32 


Optics Specifications 


Spectral Resolution 

125 nm (32 bands) or 100 nm (256 bands) 

Angular Resolution (mrad) 

.9x2 

Focal Length (microns) 

111 

Slit (microns) 

200 

Aperture (mm) 

35 


Performance Specifications 


Quantization 

12 bit 

Integration Time (ms) 

3 

Frame Rate (Hz) 

150 

Data Rate (Mbps) 

10 

Output Data Format 

AHI convertible to ENVI image format 

Display 

Flight real time desplay, SVGA Monitor 

Operating System 

PC-based, Windows NT 

Software 

Custom software for data and image acquisition, post 
mission review and processing, archival support 

Data Archive 

8mm Exabyte tape or IDE hard disk 

Data Storage 

12 GB -hot swappable 

FOV 

7 degrees cross track - pushbroom hyperspectral imager 

Guidance Sensor 

C-MIGITS II: heading, pictch, roll, lat, Ion 

Position Accuracy 

< 10 m 

NEDT 

<.1 K@300 K 

NESR 

< .02 watts/m A 2-sr-micron 

Radiometric Calibration 

NIST traceable 

Ancillary Data 

Color linescan CCD with twice the IR swath width 


Table 9. Physical and Performance Specifications of the AHI System and Sensor 


bodies, and the measurement modes. Several performance parameters of the sensor can 
be observed via plots such as noise and radiance histograms that can be displayed in real¬ 
time on the system’s computer monitor. This is illustrated in Figure 31. Additional 
details on the system and sensor physical and performance specifications are found in 
Table 9. (Williams, “Sensor Specifications”) 
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c. Data Collection Techniques and Applications 

(1) Data Collection. AHI collects data in a pushbroom fashion. 
The focal plane array has 256 rows of 256 pixels each. The 7.5 - 11.5 micron band of 
wavelengths is spread over the 256 rows so that each row represents a narrow band 
within this broad LWIR region. AHI simultaneously obtains spectra for 256 points on 
the ground. As the focal plane array is pushed along the flight line, these rows will begin 
to overlap and form a hyperspectral data cube with 256 pixels along the x-axis, a number 
of pixels along the y-axis that is determined by the duration of collection along the 
flightline, and 256 x-y image planes in the z-axis, each representing one of the narrow 
bands of wavelengths within the 7.5 - 11.5 micron region. Figure34 illustrates this 
pushbroom concept using a hyperspectral FPA and the resulting hyperspectral data cube. 


Wavelength Axis 

7.5-11.5 microns 
Range Spread Over 256 Rows 



Figure 34. Pushbroom Hyperspectral Data Collection and Resulting 

Data Cube 


Calibration is performed using the on-board adjustable temperature 
black body. High and low temperatures are chosen for the black body based upon the 
expected maximum and minimum temperatures to be imaged in the scene. A middle 
temperature for the black body is set near the average of the high and low temperatures. 
One hundred frames of data are collected at each temperature. The high and low 
temperatures are used to compute gains and offsets that are applied to data collected 
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using the middle temperature producing corresponding radiance data. Signal to noise is 
calculated over all the calibrated pixels for each wavelength so that spatial nonuniformity 
noise is included in the calculations. The black body data is also used to identify bad 
pixels in the focal plane array. A bad pixel is characterized as non-responsive, saturated, 
or “noisy” in the collection of the 100 frames of black body data. (Lucey et al., 2000, p. 
34-35) 

(2) Sensor Applications And Collection Platforms. AHI was 
originally designed for real-time mine detection. For this mission, the sensor operates in 
a pod mounted under a helicopter. Following successful demonstration of this mission, 
HIGP explored additional applications and collection platforms for AHI. AHI currently 
has three main collection platforms: helicopter, fixed wing Twin Otter or Navajo 
aircrafts, and various ground platforms. Figure 35 displays AHI operating from its three 
main collection platforms. 



Figure 35. AHI Collection Platforms: Helicopter, Ground Based, and 

Fixed Wing Aircraft 

To date, HIGP has demonstrated the following applications with 

AHI: 

• Airborne detection of land mines 

• Hyperspectral land mine phenomenology 

• Concealed target detection and phenomenology 

• Gas detection 

• Active laser hyperspectral imaging 

• Geologic mapping 

• Coastal water temperature mapping 

• Missile defense intercept test support 

• Hyperspectral basic research data collection 
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Figure 36. Single Band Images and Single Pixel Spectra Demonstrating 

AHI’s Gas Detection Capability (from Lucey et al., 2000, p. 35) 

This paper further explores AHI’s application in gas detection. An initial experiment was 
conducted in which AHI flew over a site multiple times where different known quantities 
of gas were released. Single band images and a single pixel spectrum from this 
experiment are shown in Figure 36. Figure 37 illustrates the use of a matched filter on 



Figure 37. Using Matched Filter On Single Band Images Indicating Gas 
Detection In Hashed Area (from Lucey et al., 2000, p. 36) 
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this data and confirms AHI’s ability to detect and characterize gases with broad spectral 
features. This paper extends AHI’s ability beyond the basic gas detection and 
characterization to quantifying a gas concentration. In Chapter IV, data collected over 
Hawaii’s volcanoes are used to demonstrate this extended capability through analysis of 
the volcanic SO 2 plumes. 

B. SUMMARY OF IR DETECTION 

Using COSPEC data analysis results as ground truth, Realmuto et al. successfully 
demonstrated the ability to detect and quantify SO 2 plumes using TIMS, NASA’s 
airborne multispectral LWIR sensor. In 1994, Realmuto et al.’s estimation procedure, 
based on a least squares fit between radiance profiles generated by the LWIR sensor and 
MODTRAN simulations, produced ground temperatures and SO 2 column abundances for 
single pixels in a data set. In 1997, the application of the estimation procedure was 
expanded and applied to groups of adjacent pixels, developing the ability to produce a 
two-dimensional map of an SO 2 plume. Finally in 2000, Realmuto et al., simulated the 
application of these estimation and mapping procedures to data that might be generated 
by an LWIR sensor employed on a satellite. In each of these cases, the results produced 
by the LWIR sensor reasonably compared with analysis performed by COSPEC. 

Although the LWIR sensors proved successful, Realmuto et al. did not suggest 
that they be used to replace COSPECs. It was suggested that a combination of LWIR and 
UV detection and quantification would be best. Using multiple regions of the 
electromagnetic spectrum and varying collection methods will compensate for the 
weaknesses noted in each, as well as compound their strengths. (Realmuto et al., 2000, p. 
15071) 

Realmuto et al.’s estimation and mapping procedure form a basis for data analysis 
techniques and methods to be applied to data collected by AHI in the following chapter. 
A new 2002 data set collected by both AHI and COSPEC over the Pu’u’O’o vent is used 
for comparison. 
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IV. DATA COLLECTION AND S0 2 DETECTION 


A. DATA COLLECTION 

A data collection experiment was conducted on April 18, 2002. Figure 38 
displays two pictures of the Pu’u’O’o Vent of Kilauea Volcano in Hawaii. The picture 
on the left is a visible image of the vent taken from an aerial view looking west- 
southwest. The picture on the right is an LWIR image of the vent taken by AHI mounted 
on a twin-otter plane also looking west-southwest. Both images were taken on April 18, 
2002. The AHI image is one of several that were collected for data in support of this 
paper. 



Figure 38. Visible and Infrared Image of Pu’u’O’o Vent of Kilauea 


Volcano, Hawaii, Collected on April 18, 2002 

COSPEC, FLYSPEC, and AHI collected data on the S0 2 plume of the Pu‘u‘0‘o 
vent on the morning of April 18, 2002. The original plan involved all three sensors 
collecting concurrent and corresponding data sets; however, several timing, 
communication, and weather issues prevented this from happening. 

1. Weather Conditions 

The quality of data collected via remote sensing instruments is often at the mercy 
of Mother Nature and the weather. As discussed in Chapter II, clouds and wind are two 
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of the major weather concerns in the UV analysis of SO 2 plumes. Chapter III notes that 
all aspects of the atmospheric profile between the ground and the sensor impact the IR 
analysis of SO 2 plumes. In both cases, weather conditions are an important aspect of data 
analysis and must be noted. 

COSPEC, FLYSPEC, and AHI collected data between the hours of 9:00 and 
11:00 a.m. Weather observations for these two hours were obtained from the weather 
station at the Hilo Airport, approximately 37.5 km north of Pu‘u‘0‘o and at an elevation 
of 44 feet above mean sea level. Between these hours the station reported clear skies 
with a few scattered clouds, an average surface air temperature of 26.7 degrees Celsius, a 
relative humidity of 68 percent, and a barometric pressure of 1013.5534 millibars. The 
winds were very light and variable out of the east-southeast. 

Most of these conditions were ideal with the exception of the winds. Ideal wind 
speeds and directions for data collection on the Pu’u’O’o plume are noted in Chapter II as 
northeasterly trade winds with a consistent speed of 5 m/s. These ideal wind conditions 
normally carry the plume from Pu’u’O’o down over the Chain of Craters Road. The east- 
southeast winds of April 18 carried the winds over Highway 11. The light and variable 
speeds prevented the analysis of an emission rate for the data collected, as a steady and 
consistent wind speed is one of the main variables in this calculation. A map of the 
routes covered by each of the sensors is displayed in Figure 39. 

2. COSPEC and FLYSPEC Data Collection 

Mounted on the same vehicle, COSPEC and FLYSPEC collected data on April 18 
using the ground mobile technique. As the light and variable winds made the location of 
the plume somewhat unpredictable, the sensors made an initial collection down the Chain 
of Craters Road looking for any possible signs of SO 2 . There were no indications of the 
plume along the Chain of Craters Road that morning. 

Collections were then made along a portion of Crater Rim Drive and Highway 11. 
The sensors collected data for nearly 30 minutes (9:58-10:28 a.m.) traveling at a constant 
speed of approximately 30 mph in the northeast direction along Highway 11. A four- 
minute section of the data collected in this traverse indicates a portion of plume 
containing concentrations of SO 2 as high as 747 ppm-m. This 30-minute route was 
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Figure 39. 
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Routes of Data Collection for COSPEC, FLYSPEC, and AHI on 
April 18, 2002 


retraced (10:30-10:56 a.m.) traveling approximately 30 mph in a southwest direction and 
produced another four-minute section of data indicating a section of plume with 
concentrations of SO 2 as high as 705 ppm-m. The routes covered by COSPEC and 
FLYSPEC are highlighted in black with black text indicating times and direction of travel 
in Figure 39. Additional details on the data analysis results are presented later in the 
chapter. 

3. AHI Data Collection 

AHI, mounted in a Piper Navaho plane flying between 4000 and 4500 feet above 
sea level, collected data over three main areas. Approximately 3 minutes (9:34-9:37 
a.m.) of data was collected over the chain of craters road followed by 3 minutes (9:44- 
9:47 a.m.) of data along Highway 11. Two minutes (9:52-9:54) of data were collected 
over the Pu’u’O’o Vent followed by a final six-minute collection (10:00-10:06) over 
Highway 11. The routes covered by AHI are highlighted in red with red text indicating 
times and direction of travel in Figure 2. 

During AHI’s flight, multiple files of data were also collected imaging the on¬ 
board black body for calibration purposes. Black body temperatures covering a range 
from 25 to 50 degrees Celsius were imaged. These black body images were used to 


53 






convert the digital numbers registered by AHI to values of radiance with units of 
W/m 2 /sr/|X. 


As the AHI data were calibrated, each of the runs was spatially segmented into 
smaller files to facilitate the analysis. The segments were also spectrally binned to 
reduce the number of bands from 256 to 205. Table 10 lists the four data collections, 
their associated run names, their corresponding segments, flight times, and image sizes. 
Each of the data collections will be referred to by its run name and the starting flight time 
of the appropriate segment from this point on (e.g. R1_A 093134). 


Data Collection 

Run Name 

Segement 

Starting Flight Time 
(hh:mm:ss in HST) 

Image Size 

(samples x lines x bands) 

Chain of Craters Road 

R1_A 

1 

09:31:34 

256x1499x205 



2 

09:32:26 

256x1499x205 



3 

09:34:33 

256x1499x205 



4 

09:35:25 

256x1499x205 



5 

09:36:18 

256x1499x205 



6 

09:37:18 

256x631 x 205 

First Highway 11 

R1_B 

1 

09:44:50 

256x1499x205 



2 

09:45:42 

256x1499x205 



3 

09:46:35 

256x1499x205 

Pu'u'O'o Vent 

R1_C 

1 

09:52:58 

256x1874x205 



2 

09:53:49 

256x1436x205 

Second Highway 11 

R2_B 

1 

10:00:30 

256x1499x205 



2 

10:01:21 

256x1499x205 



3 

10:02:13 

256x1499x205 



4 

10:03:12 

256x1499x205 



5 

10:04:03 

256x1499x205 



6 

10:04:55 

256x1499x205 



7 

10:05:49 

256x1499x205 


Table 10. Summary of Files Produced From AHI Data Collection 


B. S0 2 DETECTION 

The analysis of AHI’s data was approached with two basic goals in mind. The 
first goal was to demonstrate AHI’s ability to detect the presence of SO 2 . The second 
goal was to attempt to quantify the amount of SO 2 detected. This chapter addresses SO 2 
detection, while Chapter V addresses the quantification of SO 2 . 

The analysis performed to determine AHI’s ability to detect SO 2 can be broken 
down into three main areas: initial detection, application of spectral classification tools, 
and the development of spectral libraries. The data were examined in both radiance and 
emissive space in each of these areas. 
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1 . 


Initial Detection 


The first area of analysis involved an initial display of the 18 segments of AHI 
data using the decorrelation stretch method applied by Realmuto et al. to TIMS data. As 
discussed in Chapter III, the decorrelation stretch is a transform that helps remove the 
high correlation between bands of a multispectral dataset. A false color RGB image can 
then be produced to highlight areas of interest. In Realmuto et al.’s specific application 
of this technique to SO 2 detection, the color blue is assigned to the TIMS Channel 2 (8.6 
- 9.0 micron band) that corresponds to a prominent absorption feature of SO 2 . Any SO 2 
that is present in one of these false color radiance images appears in shades of yellow and 
red as the blue color representing the radiance in Channel 2 is absorbed. 


Each of the AHI radiance images was spectrally resampled reducing the 205 
bands to 6, simulating the six spectral bands of TIMS. A decorrelation stretch was 
applied to each of the simulated TIMS radiance images producing false color RGB 
images to be examined for any significant areas of yellow color. Appendix A includes 
the 18 false color RGB images created in this phase and their corresponding AHI 
radiance image. Only one image appears to have any significant yellow areas. The false 



Decorrelation Stretch - False Color RGB - Radiance 



Figure 40. False Color RGB Radiance Image Generated by the Decorrelation 

Stretch Method 


color RGB image of R1C 095349, a segment of data collected over the Pu’u’O’o Vent, 
is pictured in Figure 40 with its corresponding AHI radiance image. Several radiance 
spectra of yellow pixels are compared with radiance spectra of non-yellow pixels in 
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Figure 41 in order to verify the presence of the SO 2 absorption feature in the yellow 
pixels. The locations of the pixels used for this verification are identified in Figure 40. 
R1C095349 is the main focus for the rest of the analysis. 


SO, Spectra 



Figure 41. Radiance Spectra of Yellow Pixels Demonstrating the SO 2 

Absorption Feature 

2. Application of Spectral Classification Tools 

The second area of analysis involved applying several different spectral analysis 
tools to R1C095349. The application of a Forward Principal Component Transform, 
Spectral Angle Mapper (SAM) Method, and Matched Filter (MF) Mapping Method to the 
image is discussed in the following paragraphs. Results produced by all of these spectral 
analysis tools are presented in both radiance and emissive space. Interesting comparisons 
are presented between not only the radiance and emissive spaces, but also between the 
analysis tools. 
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a. Conversion From Radiance to Emissive Space 

In order to convert the AHI radiance image to an emissive image, the well- 
known problem of separating temperature and emissivity in thermal infrared data must be 
addressed. Kealy and Hook describe the spectral emissivity of a material as a measure of 
its ability to emit radiation compared to a blackbody. In Equation 4.1 emissivity is 
defined as the ratio of the radiance of a material that may be observed by AHI (Lahi) to 
that of a blackbody at the same temperature. Lahi can be fully described by the 
MODTRAN radiance equation described in Chapter III. Equation 4.2 is the Planck 
function that describes the radiance of a black body (L B b) in terms of the temperature 
being imaged and the wavelength used for imaging. With remotely sensed data, there are 
two unknowns in Equation 4.1, the spectral emissivity and the temperature of the 
observed material. Several methods have been developed to try and solve this 
underdetermined problem. Kealy and Hook evaluate three of the most common methods 


£(A) = 


L AM T ,) 

L bb (A,T 0 ) 


(4.1) 


e (X) = Spectral emissivity of material as a function of the wavelength 
used for imaging 

Lahi (X, T 0 ) = Radiance of material observed by AHI as a function of the 
wavelength used for imaging and the temperature of the material observed 


L B b (X, T 0 ) = Radiance of a black body as a function of the wavelength 
used for imaging and the assumed true temperature of the material being 
observed 


L bb (A,T 0 ) = 




?i 71 


f c 7 ) 


exp 


-1 






(4.2) 


Ci = First radiation constant = 3.74151 x 10~ 16 Wm 2 
C 2 = Second radiation constant = .0143879 mK 
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used and suggest that the emissivity normalization method is one of the more accurate 
methods to use particularly when working with scenes dominated by rocks and soils. 
(Kealy, P. and Hook, S., 1993, p. 1163) Realmuto successfully demonstrated the 
application of this technique to mapping the basalt flows in Hawaii. Given this 
background, the emissivity normalization method was chosen to convert the AH1 
radiance image to emissive space. 


The emissivity normalization method starts by assuming a constant value 
for emissivity. A value of .96 is commonly used and is chosen for this particular case. 
With Lahi and the emissivity defined, Equation 4.1 can be arranged to solve for the 
temperature of the material being observed as seen in Equation 4.3. The temperature of 
each pixel is calculated for each of the 205 wavelengths, and the highest temperature 
calculated is chosen as the temperature of that pixel. This temperature is then used as a 
known value in Equation 4.1 to calculate the emissivity of the pixels at each wavelength. 



Figure 42 displays the emissivity and temperature images of R1_C_095349 calculated 
using this normalized emissivity approach. 

As discussed in Chapter 111, the observed ground temperature is one of the 
important unknowns in the underdetermined problem of solving for SO 2 plume 
concentration. The importance of the temperatures in this problem warranted developing 
a means to verify the results calculated using the normalized emissivity approach. 
Several random pixels from both the radiance and emissivity images of R1C095349 
were chosen to verify the calculation. These pixels are annotated in Figure 42. Band 35, 
corresponding to a wavelength of 8.5344 p and centered on the SO 2 absorption feature, is 
used to present the radiance and emissivity images in Figure 42. Apparent instrument 
artifacts are noted along the edges of the emissivity image. 

The calculated temperature of each of the chosen pixels and an assumed 
emissivity value of .96 were reapplied to Equation (1), generating a simulated AH1 black 
body radiance profile that does not include any atmospheric or SO 2 effects. This 
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R1C 095349 - Radiance R1C095349 - Emissivity R1C095349 - Temperature 



Figure 42. R1C095349 Emissivity and Temperature Images Generated 

By the Normalized Emissivity Method and the Corresponding Radiance Image 


simulated black body profile was then plotted against the true AHI radiance profile for 
each pixel. Figure 43 displays both the radiance comparison and corresponding 
emissivity profiles generated for a 3 x 3 pixel average around the original pixels 
annotated in the images of Figure 42. The profiles representing AHI data are depicted 
with both a line and symbol whereas the simulated blackbody profiles are depicted with a 
solid line only. As discussed in Chapter III, the longer wavelengths (approximately the 




Figure 43. R1C095349 Black Body and AHI Radiance Comparison 

With Corresponding Emissivity Plot 


9.5 to 11.5 |i range) are not as effected by the atmosphere and SO 2 , so this is where a 
good fit of the two profiles is desired. Figure 43 demonstrates a reasonably good fit 
between the radiance profiles within these wavelengths. 
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The importance of understanding the effects of the method used to 
separate the temperature and emissivity was great enough to test an additional segment 
from the dataset. Figure R1_A_093718 was chosen due to the presence of water in the 
segment. As water has proven to act as a blackbody in LWIR wavelengths, the 
normalized emissivity method was applied to this segment with an assumed emissivity of 
1.0. The resulting temperatures from several pixels of both water and lava were again 
used to compare the normalized emissivity results to those generated with a simulated 
AHI black body curve. Figure 44 displays the temperature and emissivity image results 
from the normalized emissivity method. Band 35, corresponding to a wavelength of 
8.5344 p and centered on the SO 2 absorption feature, is used to present the radiance and 
emissivity images. The pixels used for the black body comparison are annotated. Figure 
45 displays the radiance comparison and emissivity profiles for the chosen lava pixels. 
The AHI data is again shown in symbols and lines while the simulated black body data is 


R1A 093718 - Radiance R1_A_093718 - Emissivity 



R1A 093718 - Temperature 


(86, 442) 

A 

(154,451) 



Figure 44. R1_A_093718 Emissivity and Temperature Images Generated 

By the Normalized Emissivity Method and the Corresponding Radiance Image 

shown in a solid line only. In this case, the comparison produced some interesting 
results. One of the lava profiles (pixel 154, 451) demonstrated a reasonably good fit with 
the simulated black body, but the other profile (pixel 86, 442) clearly did not. This is a 
result of the fact that the normalized emissivity method chooses the temperature of a 
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Lava Spectra - Run 1A 




Figure 45. R1A 093718 Black Body and AHI Radiance Comparison of Lava 

Pixels With Corresponding Emissivity Plot 

pixel to be the highest temperature calculated in the method. There appear to be several 
anomalous data points near the 11 micron wavelength that force the temperature of the 
pixel to appear greater than it is in reality. The application of this falsely inflated 
temperature forces the blackbody profile to display significantly higher values than the 
data profile. 

Figure 46 displays the radiance comparison profile for a 3 x 3 pixel 
average around the chosen water pixels. A similar effect of anomalous pixels is again 
readily seen in one of the profiles (pixel 217, 612). The anomaly of the second profile is 



Figure 46. R1A093718 Black Body and AHI Radiance Comparison of 

Water Pixels 
61 






not as easily explained. Kealy and Hook do not recommend the normalized emissivity 
method for scenes which may contain a mixture of rocks/soils and vegetation/water. 
They conducted several experiments that showed the normalized emissivity method could 
be incorrect by as much as 2.02 degrees K in a case such as this one. (Kealy, P. and 
Hook, S., 1993, p. 1163) Some of the anomalies seen may be a direct result of the 
method used. 

These comparisons provided additional insight to the effects of applying 
the normalized emissivity method to calculate the temperature and emissivity of the data 
segments. In general, it proved to be a reliable method for the R1C095349 image; 
however, it did not prove to be without fault. Keeping the possibility of anomalies in 
mind, the analysis of the data set was continued using the temperature and emissivity 
results generated by the normalized emissivity method. 

b. Application of Spectral Classification Tools 

The next area of analysis was to further develop knowledge gained from 
applying the decorrelation stretch to both the radiance and emissivity images. The goal 
was to use a spectral mapping tool to locate and map the other pixels within the images 
that displayed the SO 2 absorption feature. Two different mapping tools were applied in 
both radiance and emissive space: the Spectral Angle Mapper (SAM) and the Matched 
Filter (MF). 

A principal component analysis (PCA) was applied to both the radiance 
and emissivity images before applying either the SAM or MF. A PCA is a very common 
method of decorrelating data that requires no a priori knowledge of the scene. The 
application of a forward principal component rotation produces a user-defined number of 
new image bands called principal components (PCs). These new bands generally do not 
have any real physical meaning, but represent linear combinations of the original spectral 
bands with enhanced contrast and variance. (Stefanou, M., 1997, p.48) Eight ROI’s were 
selected using both the newly transformed PC image and the false color RGB TIMS 
image. These ROI’s included pixels that were believed to be SO 2 and those that were not 
and are shown in the radiance image in Figure 47. These ROI’s were used as endmember 
spectra for both the SAM and MF analysis. 
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The first algorithm which was then applied was the Spectral Angle 
Mapper (SAM), which operates by taking the dot product of each spectrum with a 
sequence of reference spectra. In this case, the reference spectra are defined by the ROI’s 


R1 C 095439-ROI’s 



Figure 47. Regions of Interest Used in the Application of Spectral 

Mapping Tools 


chosen from the PC images. These reference spectra and the image pixels are treated as 
vectors. The SAM compares the angle between these two vectors in a space with 
dimensionality equal to the number of image bands. If the angle between the vectors is 
small, the vectors are considered matched and of a similar class. A threshold angle can 
be chosen, and if the pixel vector does not fall within this angle in comparison to any of 
the reference vectors, the pixel will not be classified. (Research Systems Inc, 2001, p. 
511) 

The output of the SAM includes an all-encompassing color image that is 
classified according to the all of the chosen reference spectra and a number of gray-scale 
rule images that correspond to each of the individual reference spectra. The SAM, with a 
threshold angle of .10 radians, was applied to both the radiance and emissivity 
R1C095349 images employing the same eight ROI’s for each. Figure 48 displays the 
classification images for radiance and emissivity. The black strips on the edges 
correspond to the anomalous pixel regions found in the emissivity displays illustrated 
above. In much of the analysis proceeding from this point, processing was restricted to 
an 800-line subset of the run, since there was no SO 2 outside this region. The red and 
orange areas in the classification results correspond to two classes of SO 2 and will be the 
primary focus for additional analysis. The other colors in the image correspond to 
regions with no obvious SO 2 signature. There are slight differences in the SAM results 
from the two versions of the data. 
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R1C 095439 R1 C_095349 

Radiance - SAM Results Emissivity - SAM Results 



Figure 48. Initial SAM Results in Radiance and Emissive Space 

Figure 49 displays the spectral angle in gray-scale form for the red and 
orange classes. The red class images are scaled from 0 to .11 radians (white to black) and 



Figure 49. SAM Rule Images for Red and Orange S02 Classes in Radiance 

and Emissive Space 

the orange class images are scaled from 0 to .15 radians (white to black). Bright or white 
areas represent the presence of SO 2 . Slight differences do appear in the rule images for 
both classes of SO 2 . Scatter plots comparing the red and orange SO 2 rule image results in 
radiance and emissive space are displayed in Figure 50. The x and y axis are in units of 
spectral angles and represent the results for the radiance and emissive spaces 
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respectively. The smaller spectral angles represent good matches for that particular class 
of SO 2 . The color contours of the plot represent the number of occurrences in each of the 
color areas. A first order linear fit is applied to the data points to show the relationship 
between the two spaces. A correlation factor is also calculated for the data. Although 
there are slight differences in the images, the results show a good linear relationship, and 
the correlation factor calculated for both classes is relatively high. As temperature plays 
a significant role in the SO 2 signature, these differences may stem from anomalies created 



Figure 50. Scatter Plots of SAM SO 2 Rules in Radiance and Emissive 
Space for the Red and Orange SO 2 Classes, Correlation Coefficient For Both Classes = 
.82, Red Class Slope = 1.10, Intercept = .01, Orange Class Slope = .95, Intercept = .02 

by the emissivity and temperature separation discussed earlier. The differences could 
also be an artifact of the classification tool used. 

A second classifier, the Matched Filter (MF), was applied to see what sorts 
of dependence the above results had on the type of classifier being used. The technique is 
a vector projection method that uses the scene statistics, and a projection operator to 
distinguish the desired spectrum from the background. It is commonly used in the signal 
processing application. (Stefanou, M., 1997, p.87) 

In this case, the MF is used to find matches in the scene for the eight 

ROI’s employed in the SAM analysis. Each of the image pixels is compared to the 

reference ROI spectrum and a relative degree of match is assigned to each pixel. A 

perfect match is given the value of 1. The results for the red and orange classes are 

shown in Figure 51. Again, the classifier results have been scaled from dark to light, 
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Figure 51. MF Gray-Scale Images for Orange and Red Classes of SO 2 in 

Radiance and Emissive Space 

with high probability corresponding to white. The noticeable differences seen in the 
images between the radiance and emissive spaces using the SAM do not readily appear in 
the MF. The scatter plots displayed in Figure 52 support the apparent similarities 


Matched Filter - 04/18/2002 - 095349 



Figure 52. Scatter Plots of MF SO 2 Results in Radiance and Emissive Space 
for the Red and Orange SO 2 Classes, Red Class Correlation Coefficient = 0.94, 
Slope = 0.99, Intercept = 0, Orange Class Correlation Coefficient = .90, 

Slope = 0.96, Intercept = 0 
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between the classifications in radiance and emissive spaces using the MF. The x and y 
axis of the scatter plots are in units of “relative match” with 1 representing the best match 
to the orange or red class SO 2 reference spectra. The data show a strong linear 
relationship between the radiance and emissive space with a correlation factor that is 
almost 10% greater than that of the SAM results. 

A final set of scatter plots are displayed to help directly compare and 
visualize this apparent difference in the performances of the two spectral mapping tools. 
The SO 2 SAM Rule results are plotted against the SO 2 MF results for both the radiance 
and emissive spaces. Figure 53 displays the scatter plots for both SO 2 classes in radiance 



Spectral Angle (degrees) 


Radiance-04/18/2002- 

"Red Class" 

1.5“ 


095349 

S0 2 "Orange Class" 


Spectral Angle (degrees) 



Figure 53. Scatter Plot of SAM and MF Results for SO 2 Classes in 
Radiance Space, Red Class Correlation Coefficient = -0.49, Slope = -0.20, 
Intercept = 0.82, Orange Class Correlation Coefficient = -0.44, 

Slope = -0.13, Intercept = 0.54 


space. The best matches for SO 2 form the diagonal feature in the left portion of the 
scatter plot. These pixels are closest to 1 on the MF scale and have the smaller angles on 
the SAM scale. The color contours of the plot again indicate the number of occurrences 
within each of the color areas. The dense clusters on the right side of the plot are all of 
the occurrences of pixels not containing any SO 2 . A first order linear fit is applied to the 
SO 2 pixels on the left side of the plot to determine the relationship between the SAM and 
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MF results. Although the data for both classes does maintain somewhat of a linear trend, 
the correlation factors are only at -.49 and -.44. It is also interesting to note the difference 
between the red and orange classes. These results support the differences in correlation 
factors between the radiance and emissive spaces in the earlier plots. 

Figure 54 displays the scatter plot of SAM and MF results in emissive 
space for both the red and orange classes of SO 2 . A first order linear fit is applied to the 
SO 2 pixels on the left side of the plot to describe the relationship between the SAM and 
MF results. The data shows somewhat of a linear trend, although the correlation factor is 
lower than it was in radiance space at -.37 and -.39. It is again interesting to note the 
differences between the two classes. The red class consistently produces higher 
correlation factors in each of the scatter plots. The differences seen in the spectral 
mapping of the red and orange classes encouraged additional analysis of the SO 2 class 
results. 


Emissivity - 04/18/2002 - 095349 


S0 2 "Red Class" 



MF = 0.71 -I- -0.18 *SA 

Correlation = -0.37 
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MF = 0.52 + -0.13 *SA 

Correlation = -0.39 
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Figure 54. Scatter Plot of SAM and MF Results for SO 2 Classes in 
Emissive Space, Red Class Correlation Coefficient = -0.37, Slope = -0.18 
Orange Class Correlation Coefficient = -0.39, Slope = -0.13, Intercept = 0.52 

3. Development of Spectral Libraries. 

In addition to using the spectral classification tools, the red and orange SO 2 
classes are examined in Environment for Visualizing Images’ (ENVI’s) n-Dimenaional 
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(n-D) Visualizer with hopes of improving the understanding of the differences made 
apparent in the SAM and MF analysis. The n-D Visualizer allows a user to create scatter 
plots of selected data in n-dimensions where “n” is the number of bands. The user can 
interactively rotate selected data in n-Dimensions gaining new insight to the composition 
of the data, with the opportunity of refining selected classes or ROI’s. 

The initial red and orange SO 2 classes were examined in both radiance and 
emissive space. Using the PC band radiance and emissivity images, the red and orange 
SAM SO 2 ROI data were exported to the n-D Visualizer. Six spectrally distinct classes 
of SO 2 were developed. Each class was assigned a different color and had different 
numbers of pixels that belonged to them. Table 11 is a summary of the different classes 
and states how many pixels were assigned to each through the analysis process. 

Figure 55 is a 2-D scatter plot using PC bands 1 and 2 in radiance space. The x- 
axis of the scatter plot, corresponding to PC band 1, represents average brightness with 
brightest pixels in the ROI on the right. The upper left hand comer shows a cluster of 
blue pixels representing the coolest and darkest class of SO 2 pixels in the scene. The 
green and cyan pixels on the right end of the scatter plot represent the warmest and 
brightest classes of SO 2 pixels in the scene. 


S0 2 Spectra 



Figure 55. 2-D Scatter Plot of S02 Classes in Radiance Space Using 

PC Bands 1 and 2 and the Corresponding Radiance Spectral Library of the SO 2 Classes 
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A spectral library was generated from these six classes of SO 2 . The spectral 
library helps illustrate the differences between the classes in numeric form and maintain 
records for future analysis. The profiles for the SO 2 radiance spectral library are 
displayed in Figure 55. The profiles reinforce the fact that there is clear separation 
between each of the classes. The brightest or warmest pixels display radiance values as 
high as 25 to 33 W/m 2 /sr/|0, while the darkest or coolest pixels display radiance values 
between 11 and 13 W/m 2 /sr/|0,. 

Figure 56 is a 5-D scatter plot using the first five PC bands in emissivity 
space. The distinction between the classes is more difficult to make using the n-D 
Visualizer in emissivity space; however, the spectral library profiles, also displayed in 
Figure 56, again demonstrate a clear separation of the SO 2 classes. It is interesting to 
note where the classes of pixels fall in the emissivity space. The warmest and coolest 
classes in radiance space also fall at either end of the emissivity plot. However, the 
middle region displays information worth noting. In radiance space, the green class of 


S0 2 Spectra 



Wavelenqth(u) 

Figure 56. 5-D Scatter Plot of S02 Classes in Emissivity Space Using the 

First Five PC Bands and the Corresponding Emissivity Spectral Library of the SO 2 

Classes 



pixels appear to be a warmer subset of the red class. In emissivity space, the green class 
appears to have emissive characteristics similar to the orange class. The coral class 
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appears to be a subset of the orange class in radiance space. While in emissive space, it 
takes on characteristics similar to the red class. 


In both emissive and radiance space, the prominent difference noted between the 
classes is a y-axis offset from one another. This offset is partially explained in 
temperature differences. The six ROI classes were applied to the R1C095349 
temperature image and statistics were computed for each class. These statistics are 
compiled in Table 11. Although there is some overlap in temperature ranges between the 


Class Color 

Class Size (points) 

Min Temp (K) 

Max Temp (K) 

Mean Temp (K) 

Stdev 

Green 

64 

384.41 

425.41 

398.84 

12.01 

Cyan 

649 

357.12 

425.80 

380.02 

13.71 

Red 

7226 

318.81 

398.07 

350.08 

14.62 

Orange 

4083 

319.72 

377.34 

345.09 

11.64 

Orchid 

1911 

314.28 

373.10 

329.07 

11.30 

Blue 

32 

311.80 

319.11 

315.19 

1.99 

Table ] 

1. SO 2 Class Temperature Statistics 


classes, each class has a distinct mean. The orange and red classes show the smallest 
difference between classes with a 5-degree separation between their temperature means. 
The red class displays the largest range of temperatures and largest standard deviation, 



Figure 57. R1C095349 Radiance Image With Final SO 2 

Classes Annotated 
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while the blue class displays the smallest temperature range and smallest standard 
deviation. These results further emphasize the importance of the temperature and 
emissivity determination. The temperature has a significant impact on the spectra 
observed when SO 2 is present. Finally Figure 57 displays the six SO 2 classes on a 
radiance image of R1C095349. It should be noted that the coral class described in the 
spectral library figures is annotated as orchid in Figure 57 and in all future analysis. 
Figure 57 is the summary figure for the first goal of the AHI data analysis displaying the 
results of the SO 2 detection. These six classes are used as a starting point for the second 
goal of attempting to quantify the amount of SO 2 detected. 

C. SUMMARY OF S0 2 DETECTION 

In summary, the three areas of SO 2 detection included: initial detection, applying 
spectral classification analysis tools, and developing a spectral library of SO 2 profiles. 
SO 2 was initially detected by applying a decorrelation stretch to AHI data that was 
resampled to create a simulated TIMS image. This produced a false color RGB image in 
which any presence of SO 2 was identified by yellow colored pixels. The decorrelation 
stretch image was used as a visual aid to select several ROI’s in the image chosen for 
analysis. Some ROI’s included SO 2 while others did not. Applying the SAM and MF to 
the ROI’s created two distinct classes of SO 2 represented in colors of orange and red. 
The red and orange classes were then further analyzed in ENVI’s n-D Visualizer 
resulting in six final classes of SO 2 that are used as a starting point for the second goal of 
attempting to quantify the amount of SO 2 detected. 
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V. QUANTIFICATION OF S0 2 DETECTIONS 


As stated in Chapter IV, the analysis of AHI data was approached with two goals 
in mind. The first goal, addressed in Chapter IV, detected an S02 plume in one of the 
images collected and classified the plume into six spectrally distinct classes of S02. The 
second goal, to be addressed in this chapter, is to quantify the amount of S02 in these 
S02 classes. While detection was performed in both radiance and emissive space, 
radiance space is the primary focus for quantification. 

As discussed in Chapter III, the process used by Realmuto et al. to determine the 
amount of S0 2 in the plume is based on an underdetermined radiance equation. The 
equation is given in Chapter III and has two unknown variables of temperature and S0 2 
concentration. The method used to solve this equation for the AHI data is similar to the 
method used by Realmuto in that it is based on a comparison of the sensor generated 
radiance profdes and radiance profiles generated by MODTRAN. Details of the method 
are discussed in the following paragraphs. 

A. MODTRAN 

1. Inputs 

The basic inputs for the MODTRAN code consist of four main “cards” or lines of 
FORTRAN code. Each card defines a different portion of the radiance equation. The 
first card addresses the atmospheric conditions. MODTRAN has six model atmospheric 
profiles that can be used for inputs. They are specified according to latitudinal sections 
of the globe and season of the year. Users can also define their own atmospheric profile 
for a specific situation. A combination of these two options is used in the AHI data 
analysis. The first card also contains an input for the ground surface temperature that has 
significant impact on the analysis problem. The temperatures determined by the 
normalized emissivity method are the temperatures initially applied to this variable for 
analysis. 

The second card allows the user to define an atmospheric profile. A number of 
layers for the user-defined profile are chosen, and each layer is described according to 
temperature, pressure, relative humidity, and atmospheric constituents, which include 
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SO 2 . For the AHI data analysis, four layers were developed. The initial layer is 
developed from the ground to the altitude of the top of the plume. The temperature, 
pressure, and relative humidity of this layer are taken from the weather reports from the 
Hilo weather station during the hours of collection. This is also the layer where the SO 2 
concentration of the plume is defined. Two important assumptions relating to the SO 2 
plume are made in the development of this layer. First, an assumption is made regarding 
the altitude and the thickness of the plume. As R1C095349 is an image collected 
directly over the vent, it is assumed that the altitude of the plume is at the surface of the 
vent or the very bottom layer in the atmospheric profile. As no direct measurements on 
the thickness of the plume were collected, an estimate of 150 m is used. Second, it is 
assumed that the temperature of the plume is in thermal equilibrium with the particular 
layer of the atmosphere in which it is described. The remaining three layers of the model 
are divided between the top of the plume and the altitude of the sensor. The standard 
atmospheric model is used to define the temperature, pressure, relative humidity, and 
atmospheric constituents of these layers. An example set of cards with all of the input 
parameters used in the AHI data analysis is found in Appendix B. 

2. Outputs 

MODTRAN generates several output files detailing the model inputs and the 
solution to the radiance equation from the ground to the sensor. The output file 
describing the radiation observed at the model sensor is used for the comparison to the 
AHI data. This output file consists of the range of operating wavelengths for the model 
sensor and the corresponding radiance seen at each of the wavelengths. The first output 
generated with the model used for the AHI analysis is displayed in Figure 58. It is the 
default atmospheric radiance profile with a ground surface temperature of 288 K and no 
SO 2 applied. 

3. MODTRAN Model Examples 

Prior to starting any comparisons with AHI data, several test models were run 
modifying the general input parameters to be used in analysis. Many of the variances and 
sensitivities tested are described in Chapter III such as plume altitude, plume thickness, 
and model atmospheres. There were no significant differences in the sensitivity results of 
the model to report. 
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Standard Atmosphere Radiance Profile -288 K 



Wavelength (|x) 

Figure 58. MODTRAN Standard Atmospheric Profile 

During analysis, the only two input parameters of the MODTRAN code that are 
modified are the ground surface temperature and the SO 2 concentration. Several test runs 
were conducted varying the SO 2 concentrations for two ground surface temperatures (350 
K and 380 K) chosen from the mean temperature statistics of the six SO 2 classes 
identified in Table 11. SO 2 concentrations from 0 to 1600 ppm in intervals of 400 ppm 
were applied to both temperatures. The results for both temperatures are displayed in 
Figure 59. The resolution of the general MODTRAN model is much higher than the 
resolution of AHI. The individual data points are included in the figures; however, a 



Figure 59. MODTRAN Radiance and Relative Radiance Model Profiles 
for Multiple Concentrations of SO 2 at 350 K and 380 K 
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running average, displayed as a solid line, is applied to each of the profdes in order to 
simplify the display and discussion of the data. In addition to the general radiance 
profile, dividing each of the profiles containing SO 2 by the SO 2 free profile also 
generated relative radiance profiles. It is interesting to compare the relative radiance 
profiles for the two ground temperatures. The higher temperature of 380 K produces 
relative radiances of 5 to 10 percent lower than the temperature of 350 K with equivalent 
concentrations of SO 2 . In other words, with higher ground temperatures, less SO 2 is 
required to produce the same amount of absorption seen at lower ground temperatures. 
This is a direct demonstration that reemphasizes the importance of the ground 
temperature and emissivity separation conducted earlier in the analysis process. This 
exercise also improves the understanding of the model and its sensitivities to the varying 
input parameters, and establishes a reference point for the beginning analysis. 

B. AHI DATA CHOSEN FOR ANALYSIS 

Several pixel block areas from each of the six classes of SO 2 were chosen for 
analysis. Initially, individual pixels were chosen based on ground temperature and 
location. The determined ground temperatures of the chosen pixels are within one 
standard deviation of the mean for the class, according to the statistics calculated in Table 
11. They are also located in an area where they are surrounded by pixels of the same 
class. This provided the opportunity to form blocks of same class pixels to be used in 
generating a pixel averaged radiance profile. These pixel averaged radiance profiles are 
used in the analysis comparison with the MODTRAN radiance profiles. Table 12 lists 
the pixel blocks chosen for analysis, along with the individual pixels’ temperature 
determined through the normailized emissivity temperature and emissivity separation 
analysis method in Chapter IV. The individual pixels’ temperatures provided a starting 
point for the MODTRAN analysis. Figure 60 displays a spatially subsetted image of 
R1C095349 focusing on the six classes of SO 2 in the plume. The comer x-y 
coordinates of the chosen pixel blocks are annotated. 

A MODTRAN profile was generated using the pixel ground temperatures and an 

estimated SO 2 concentration. The MODTRAN output was plotted against the AHI pixel 

block profiles. The profiles were initially examined for fit in the longer wavelength 

region (9.5 - 11.5 p). This region is unaffected by SO 2 and can be used to examine the 
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profile in terms of temperature. Adjustments to the ground surface temperature input 
were made if an ideal fit was not observed. After a proper fit was established in the 


Class Color 

Class Size (points) 

Pixel Block Size 

Pixel Block Location 

Class Mean Temp 

Emiss Norm Temp (K) 

Red 

7226 

5x5 

87, 534 

350.08 

350.98 

Red 

7226 

5x5 

161, 603 

350.08 

350.56 

Red 

7226 

5x5 

52, 635 

350.08 

350.33 

Red 

7226 

5x5 

51, 660 

350.08 

348.40 

Red 

7226 

5x5 

88, 638 

350.08 

346.49 

Red 

7226 

5x5 

154, 637 

350.08 

337.72 

Red 

7226 

5x5 

39, 528 

350.08 

345.00 

Red 

7226 

5x5 

115, 534 

350.08 

367.61 

Red 

7226 

5x5 

166, 528 

350.08 

364.59 

Orange 

4083 

5x5 

60, 503 

345.09 

350.47 

Orange 

4083 

5x5 

141, 513 

345.09 

344.38 

Orange 

4083 

5x5 

98, 513 

345.09 

344.95 

Orange 

4083 

5x5 

66, 648 

345.09 

344.72 

Orange 

4083 

5x5 

103, 646 

345.09 

352.89 

Orange 

4083 

5x5 

81, 646 

345.09 

354.34 

Orange 

4083 

5x5 

35, 646 

345.09 

356.48 

Orchid 

1911 

6x3 

140, 664 

329.07 

329.77 

Orchid 

1911 

6x3 

210, 687 

329.07 

329.73 

Orchid 

1911 

6x3 

84, 671 

329.07 

318.47 

Cyan 

649 

6x3 

38,517 

380.02 

380.27 

Cyan 

649 

6x3 

145, 511 

380.02 

382.28 

Cyan 

649 

6x3 

137, 521 

380.02 

388.63 

Green 

64 

4x1 

83, 522 

398.84 

403.11 

Green 

64 

8x1 

212, 503 

398.84 

396.81 

Blue 

32 

2x1 

184, 653 

315.19 

315.04 

Blue 

32 

2x1 

199, 649 

315.19 

313.34 

Blue 

32 

1 x 1 

99, 660 

315.19 

317.74 


Table 12. 

AHI Pixel 

Data Chosen for Analysis 



Figure 60. X-Y Coordinates of Pixel Blocks in the Plume Chosen 

for Analysis 
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longer wavelengths, the SO 2 concentration was adjusted to develop the best fit to the SO 2 
signature seen in the AHI data. Comparison plots were generated for each of the pixel 
locations indicated in Table 12. 

C. AHI DATA ANALYSIS RESULTS 

As shown in Figure 60, there are two main sections of SO 2 . They are referred to 
as the upper and lower plume in the analysis discussion. The analysis results provide 
interesting comparisons both within and between the classes. The results are initially 
presented by class. A table highlighting maximums and minimums is presented in the 
individual sections along with a portion of the comparison plots that were generated. The 
MODTRAN model has a much higher spectral resolution than the AHI data so a 20-point 
running average is applied to all of the MODTRAN results. This is believed to contribute 
to a small leftward offset observed between the two datasets. A summary table and 
figure at the end of this section provides a combined numeric overview of all of the class 
results. A full collection of all of the plot results is found in Appendix C. 

1. SO 2 Class Results 
a. Red Class 

As the red class is the largest of the six classes, nine blocks of pixels were 
chosen for analysis. Each of the blocks is a 5 x 5 average of pixels. The red class is 


RED CLASS S02: T = 363 K S02 = 350 ppm X:166 Y:528 [5x5 AVG] RED CLASS S02: T = 335 K S02 = 550 ppm X:154 Y:637 [5x5 AVG] 




Figure 61. Upper, Lower, and Mid Plume Section AHI - MODTRAN 
Comparison Plots for the Red Class, Black = AHI data, Red = Modtran Model With 
20 Point Running Average Applied 
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found in both the upper and lower sections of the plume. It also comprises a majority of 
a small mid section of the plume. The nine pixel blocks chosen for analysis are dispersed 
throughout these three sections of plume. Plot results from three of the blocks (one from 
each section of the plume) are displayed in Figure 61. The plots demonstrated a 
reasonably good fit between the MODTRAN model and AHI data. Table 13 lists the 
maximum and minimum ground temperatures found in the red class and the sections of 
the plume in which they are found. All of the red class temperatures analyzed fall within 
one standard deviation 


Red Class S02 

deg K or ppm 

Plume Section 

Max Temp 

363 K 

Upper 

Min Temp 

335 K 

Lower 

Max S02 

725 ppm 

Lower 

Min S02 

325 ppm 

Lower 


Table 13. Maximum and Minimum Temperature and SO 2 Results for 

the Red Class 

of the class mean according to the statistics computed earlier. Although a large range of 
temperatures is created by the maximum and minimum results, a majority of the 
temperatures actually fell between 345 and 355 K. The ground temperatures in the upper 
plume appeared slightly higher than those in the lower section of plume. There did not 
appear to be a particular trend with the SO 2 in either section of the plume. The SO 2 
found in the red class cover ranges spanning 300 to 400 ppm in the upper and lower 
sections of the plume respectively. Comparing all of the class results, the red class 
displays the lowest values in SO 2 concentration. 
b. Orange Class 

The orange class is roughly half the size of the red class and has a large 
presence in both the upper and lower sections of the plume. Seven pixel blocks, each 5 x 
5 pixel averages, were plotted against the MODTRAN model. One plot from each 
section of the plume is displayed in Figure 62. The plots show reasonably good 
agreement with the MODTRAN model. Table 14 gives the maximum and minimum 
ground temperature and SO 2 results for the class. All of the temperatures analyzed fall 
within one standard deviation of the mean temperature calculated in earlier statistics. 
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Although smaller than the red class, the orange class is at least twice the size of the other 
classes, and it displayed a much smaller 

ORANGE CLASS S02: T = 343 K S02 = 1225 ppm X:141 Y:513 [5x5 AVG] ORANGE CLASS S02: T = 350 K S02 = 950 ppm X:35 Y:646 [5x5 AVG] 




Figure 62. Upper and Lower Plume Section AHI - MODTRAN Comparison 
Plots for the Orange Class, Black = AHI Data, Orange = Modtran Model Data with 
20 Point Running Average Applied 


variance in ground temperatures than the red class. Temperatures appear slightly lower 
in the upper section of plume, while the SO 2 results are significantly lower in the lower 
portion of the plume. Overall the SO 2 results also displayed a smaller variance than the 
red class with ranges spanning only 100 to 200 ppm in the lower and upper plume 
sections respectively. 


Table 14. 


Orange Class S02 

deg K or ppm 

Plume Section 

Max Temp 

350 K 

Lower 

Min Temp 

343 K 

Upper 

Max S02 

1400 ppm 

Upper 

Min S02 

850 ppm 

Lower 


Maximum and Minimum Temperature and SO 2 Results for the 
Orange Class 


c. Orchid Class 

As noted in the initial development of the classes, the orchid class appears 


ORCHID CLASS S02: T = 344 K S02 = 1100 ppm X:140 Y:664 [6x3 AVG] 


ORCHID CLASS S02: T = 323 K S02 = 950 ppm X:210 Y:687 [6x3 AVG] 




Figure 63. Lower Plume Section AHI - MODTRAN Comparison 

Plots for the Orchid Class 
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to be a cooler or darker subset of the orange class. The orchid class is roughly half the 
size of the orange class and is only found in the lower plume section. Three pixel block 
averages, each 6x3 pixel averages, were compared with the MODTRAN model and 
show a reasonably good fit. Two of the comparison plots are displayed in Figure 63. 
The orchid class shows a large ground temperature range in the results just as the red 
class did. The maximum and minimum temperature and SO 2 results are listed in Table 
15. Although all of the ground temperature results fall within the maximum and 
minimum temperature range calculated in the statistics, the maximum temperature result 
is slightly outside of one standard deviation of the mean calculated for the class. Two of 
the three orchid plots display lower ground temperatures than both the red and orange 
classes within the lower section of the plume. All of the plots display higher SO 2 results 


Table 15. 


Orchid Class S02 

deg K or ppm 

Plume Section 

Max Temp 

344 K 

Lower 

Min Temp 

317 K 

Lower 

Max S02 

1275 ppm 

Lower 

Min S02 

950 ppm 

Lower 


Maximum and IV 


inimum Temperature and SO 2 Results 


for the Orchid Class 


than the red and orange classes. The SO 2 results are closer in value to the orange class 
results in the upper plume section. This similarity reflects the fact that it initially evolved 
as a subset of the orange class. 

d. Cyan Class 

The cyan class initially evolved as a hotter or brighter subset of the orange 
class. The cyan class is a third of the size of the orchid class and is only located in the 

CYAN CLASS S02: T = 355 K S02 = 1575 ppm X:38 Y:517 [6x3 AVG] CYAN CLASS SQ2: T = 371 K SQ2 = 1150 ppm X:137 Y:521 [6x3 AVG] 




Figure 64. Upper Plume Section AHI - MODTRAN Comparison 

Plots for the Cyan Class 
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upper plume section. Three pixel block averages, each 6x3 pixel averages, are 
compared with the MODTRAN model and show a reasonably good fit. Two of the 
comparison plots are displayed in Figure 64. The cyan class has a substantially higher 
mean ground temperature than the other classes. However the range of ground 
temperatures in the results is moderate in comparison with the other classes. The 
maximum and minimum temperature and SO 2 results are displayed in Table 16. The 
temperatures of the chosen pixel blocks all fall at least 10 degrees below the calculated 
mean temperature but are still higher than those found in the orange class. The SO 2 
results display the highest values seen in all of the classes and are again on the order of 
those seen in the orange class reinforcing the fact that it evolved as a subset of the orange 
class. 


Table 16. 


Cyan Class S02 

deg K or ppm 

Plume Section 

Max Temp 

371 K 

Upper 

Min Temp 

355 K 

Upper 

Max S02 

1600 ppm 

Upper 

Min S02 

1150 ppm 

Upper 


Maximum and IV 


inimum Temperature and SO 2 Results 


for the Cyan Class 


e. Green and Blue Classes 


The green and blue classes are very small classes in comparison to the 
others and hence are presented together. The green class is found mainly in the upper 
plume section, while the blue class is found mainly in the lower plume section. They 
both appear randomly in their respective plume sections in groups of pixels numbering as 
large as eight and as small as one. 


GREEN CLASS S02: T = 388 K S02 = 550 ppm X:212 Y:503 [8x1 AVG] GREEN CLASS S02: T = 392 K S02 = 575 ppm X:83 Y:522 [4x1 AVG] 



Figure 65. Upper Plume Section AHI - MODTRAN Comparison Plots 
for the Green Class 
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The green class initially evolved as a subsection of the hottest and 
brightest pixels in the red class. Only two pixel block comparisons are made with the 
MODTRAN model and both demonstrate reasonably good fits. These plots are displayed 
in Figure 65. The temperatures are within one standard deviation of the mean 
temperature calculated in the statistics. The green class displays the highest ground 
temperature values seen in all of the classes. The SO 2 results are closer to those of the 
red class reinforcing its evolution from that class. The maximum and minimum ground 
temperature and SO 2 results for the green class are given in Table 17. 


Green Class S02 

deg K or ppm 

Plume Section 

Max Temp 

392 K 

Upper 

Min Temp 

388 K 

Upper 

Max S02 

575 ppm 

Upper 

Min S02 

550 ppm 

Upper 


Table 17. Maximum and Minimum Temperature and SO 2 Results 

for the Green Class 

The blue class is the smallest of all classes and originally evolved as the 
very coolest or darkest pixels in the scene. This is verified with the comparison plots. 
Three pixel blocks are plotted against the MODTRAN model and two of them are 
displayed in Figure 66. The comparison shows a reasonably good fit between the model 


BLUE CLASS S02: T - 316 K S02 - 800 ppm X:99 Y:660 [1x1 AVG] BLUE CLASS S02: T = 311 K S02 = 725 K X:199 Y:649 [2x1 AVG] 



Figure 66. Lower Plume Section AHI - MODTRAN Comparison 

Plots for the Blue Class 


and the data. The temperature range calculated in the statistics is very small and all of the 
results are within this range, but only one of them is within one standard deviation of the 
mean. The ground temperature results are the lowest of all the classes. The SO 2 results 
fall roughly between the red and orange classes. Although not quite as high, they appear 
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closer in value to those of the orange class. The maximum and minimum ground 
temperature and SO 2 results for the blue class are given in Table 18. 


Blue Class S02 

deg K or ppm 

Plume Section 

Max Temp 

316 K 

Lower 

Min Temp 

311 K 

Lower 

Max S02 

800 ppm 

Lower 

Min S02 

725 ppm 

Lower 


Table 18. 


Maximum and V 


inimum Temperature and SO 2 Results 


for the Blue Class 


f. Class Summary 

Figure 67 is a summary of all of the class results. It is the same picture 
shown earlier with the X-Y coordinates of the pixel blocks, but in this case, the 
coordinates are the ground surface temperature and SO 2 concentration inputs that 
produced the best fitting MODTRAN model to the data profile at each of the pixel block 
points. These results are also summarized in a table format in Table 19. 



Figure 67. Summary of MODTRAN Ground Surface Temperatures and 
SO 2 Concentrations Producing the Best Fit to Corresponding AHI Data 


As all of the results are visualized as a whole, it appears that the upper 

plume section has warmer ground surface temperatures and has higher SO 2 values than 

the lower section. The red class appears to maintain relatively consistent ground 

temperatures and SO 2 values in all sections of the plume; however this is not the case 

with the orange class. It is interesting to note that while the orange class maintains very 
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consistent ground temperatures, there appears to be a large difference in the SO 2 values 
between the upper and lower plume sections. 


Several different theories might explain the wide range in ground 
temperatures and SO 2 values seen in the data. These values may largely be an artifact of 
directly imaging the vent source. It is apparent from the visible image in the beginning of 
the chapter that the plume thickness and density is quite variable at the source. The fact 
that there was little to no wind on the day of collection increases this effect. Realmuto et 
al. saw some of this phenomenon in their studies and labeled it “puffing” within the 


Class 

Class Size 

Pixel Block 

Pixel Block 

Plume 

Class Mean 

Emiss Norm 

Modtran 

Modtran 

Color 

(points) 

Size 

Location 

Section 

Temp (K) 

Temp (K) 

Temp (K) 

S02 (ppmv) 

Red 

7226 

5x5 

87, 534 

Upper 

350.08 

350.98 

352 

350 

Red 

7226 

5x5 

161,603 

Mid 

350.08 

350.56 

345 

500 

Red 

7226 

5x5 

52,635 

Lower 

350.08 

350.33 

355 

475 

Red 

7226 

5x5 

51,660 

Lower 

350.08 

348.40 

345 

725 

Red 

7226 

5x5 

88,638 

Lower 

350.08 

346.49 

350 

325 

Red 

7226 

5x5 

154, 637 

Lower 

350.08 

337.72 

335 

550 

Red 

7226 

5x5 

39, 528 

Upper 

350.08 

345.00 

346 

650 

Red 

7226 

5x5 

115, 534 

Upper 

350.08 

367.61 

361 

450 

Red 

7226 

5x5 

166, 528 

Upper 

350.08 

364.59 

363 

350 

Orange 

4083 

5x5 

60,503 

Upper 

345.09 

350.47 

351 

625 

Orange 

4083 

5x5 

141,513 

Upper 

345.09 

344.38 

343 

1225 [ 

Orange 

4083 

5x5 

98, 513 

Upper 

345.09 

344.95 

343 

1400 

Orange 

4083 

5x5 

66,648 

Lower 

345.09 

344.72 

343 

975 

Orange 

4083 

5x5 

103, 646 

Lower 

345.09 

352.89 

346 

875 

Orange 

4083 

5x5 

81,646 

Lower 

345.09 

354.34 

349 

850 | 

Orange 

4083 

5x5 

35, 646 

Lower 

345.09 

356.48 

350 

950 | 

Orchid 

1911 

6x3 

140, 664 

Lower 

329.07 

329.77 

344 

1100 

Orchid 

1911 

6x3 

210, 687 

Lower 

329.07 

329.73 

323 

950 

Orchid 

1911 

6x3 

84, 671 

Lower 

329.07 

318.47 

317 

1275 

Cyan 

649 

6x3 

38, 517 

Upper 

380.02 

380.27 

355 

1575 

Cyan 

649 

6x3 

145, 511 

Upper 

380.02 

382.28 

362 

1600 


649 

6x3 

137, 521 

Upper 

380.02 

388.63 

371 

1150 

Green 

64 

4x1 

83, 522 

Upper 

398.84 

403.11 

392 

575 

Green 

64 

8x1 

212,503 

Upper 

398.84 

396.81 

388 

550 

Blue 

32 

2x1 

184, 653 

Lower 

315.19 

315.04 

311 

800 

Blue 

32 

2x1 

199, 649 

Lower 

315.19 

313.34 

311 

725 

Blue 

32 

1 xl 

99,660 

Lower 

315.19 

317.74 

316 

800 

Table 19. Summary of MODTRA 

N Grounc 

1 Surface 1 

temperatures and SO 2 


Concentrations Producing the Best Fit to Corresponding AHI Data 


plume. It is an artifact of being able to image a large area of the plume instantaneously. 
The large variation in ground temperatures could be attributed to large variations in 
surface composition. Without any additional data over the source or ground truth to 
compare, it is difficult to come to any additional conclusions. 

2. COSPEC / FLYSPEC Results 

Without any data collected in similar areas at similar times using similar methods, 
it is difficult to compare the IR and UV results. The absence of proper wind conditions 
denied the ability to generate any emission rates that would be worthy of documentation. 
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While no direct comparisons are drawn between the IR and UV data, the COSPEC and 
FLYSPEC results from the data collections along Hwy 11 (described in an earlier section 
of the chapter) are presented in Figures 68 and 69. It is worthy of noting a comparison 
between these two instruments for a second time. Figure 68 shows two plots of the data 
collected as COSPEC and FLYSPEC simultaneously traveled Highway 11 in a northeast 
direction. The plots are of SO 2 path concentration versus time. The figure on the left 
displays the data collected for the entire traverse. Some anomalies are observed in the 


Run3 - 4/18/02 



Seconds 


Run3 - 4/18/02 



Seconds 


Figure 68. Plots of Data Collected Simultaneously by COSPEC and 

FLYSPEC Traveling Northeast on Highway 11 from 0958 to 1028 on April 18, 2002 


COSPEC data due to saturation of the sensor when passing under trees. The FLYSPEC 
data is the original data with a 3-point smooth applied. The image on the right is an 
expanded version of the same data shown on the left between 1200 and 1700 seconds. 
Figure 69 is a plot of data collected as COSPEC and FLYSPEC retraced their path along 
Highway 11 traveling in the northwest direction. 


Run 4-4/18/02 



Figure 69. Plot of Data Collected Simultaneously by COSPEC and 
FLYSPEC Traveling Southwest on Highway 11 from 1030 to 1056 on April 18, 2002 


As first demonstrated in Chapter II, FLYSPEC again displays a notable ability to 
produce results similar to COSPEC. FLYSPEC’s results in Figures 68 and 69 were 


86 















generated on the move as data was being collected. Instead of watching data unfold on a 
mechanical plotter as done in collecting data with COSPEC, the collection output is 
viewed on the screen of a laptop. This comparison reemphasizes that the same product 
being generated over the years with COSPEC can now be generated equally well with the 
upgraded technology of FLYSPEC without losing the ability to maintain the consistency 
of any long running datasets. 

D. SUMMARY 

The detection and quantification of SO 2 plumes in the LWIR portion of the 
electromagnetic spectrum looking down from the air or space is not a trivial problem. 
The analysis is greatly affected by two very important underdetermined problems relating 
to the received radiation at the sensor. The first underdetermined problem is the 
separation of ground surface temperature and emissivity that generate the radiation seen 
at each pixel in the scene. The second underdetermined problem is finding the 
combination of the unknown ground temperature and SO 2 concentration that produces the 
radiation seen at each pixel in the scene. To complicate matters, the results of one 
underdetermined problem are used to solve the other. There are endless possible 
combinations of ground temperature and SO 2 concentrations that produce similar 
radiance profiles. As difficult as the problem appears, the methods described in this 
chapter have been developed to solve it. 
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VI. SUMMARY AND CONCLUSIONS 


This thesis examines the detection and quantification of volcanic SO 2 plumes via 
the study of UV and LWIR radiation. Discussion began with the work done in the UV 
region of the electromagnetic spectrum; as ground based UV sensors have a 30-year 
history of detection and quantification of SO 2 plumes. Work done with the University of 
Hawaii’s FLYSPEC was conducted in the traditional manner, and high sensitivity to SO 2 
was found, with good agreement between FLYSPEC and the COSPEC sensors. 

Analysis of SO 2 data collected by the UH sensor, AHI, was presented next. AHI 
was flown over the area surrounding the Pu’u’O’o Vent of Kilauea Volcano in Hawaii to 
provide the data for analysis. COSPEC and FLYSPEC also participated in the 
experiment to provide ground truth. 

The analysis is broken into two sections: the detection of SO 2 and the 
quantification of the SO 2 detected. The detection problem was approached using two 
main spectral classification and mapping tools: the SAM and the MF. Both tools were 
applied to data in both radiance and emissive space. Similar results were obtained. The 
SAM results were used to partition the data in the quantification analysis which followed, 
thus exploiting the redundant character of hyperspectral imagery. 

The emissivity normalization method was applied to separate temperature and 
emissivity. The emissivity and temperature results generated using this method are then 
used to help solve the second underdetermined radiance equation in which the two 
unknown variables are ground temperature and SO 2 concentration. 

MODTRAN was used to help solve this second underdetermined radiance 
equation by modeling the atmosphere between the ground and the sensor. Specific 
variables describing the SO 2 plume such as plume altitude, plume thickness, and 
meteorological conditions on the day of collection were used to develop a number of 
atmospheric layers describing the path traveled by radiation in the model. Specific inputs 
for the sensor such as operating altitude and operating wavelengths were also required. 
The two unspecified variables in the model were ground temperature and SO 2 
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concentration. The model was used to generate a spectral library of radiance profiles 
containing numerous variations on combinations of different ground temperatures and 
SO 2 concentrations. These MODTRAN profiles were then plotted against the AHI 
radiance profiles to determine a best fit. 

AHI reported significantly higher SO 2 concentration path lengths than the UV 
sensors, about 100 times higher. Unfortunately, the results of AHI’s IR data and 
FLYSPEC’s UV data are not directly compared due to several location, weather, and 
timing factors. There is a significant difference in the physical locations of the IR and 
UV data collections that produced SO 2 signatures. AHI detected SO 2 directly over the 
vent while FLYSPEC detected SO 2 several kilometers away from the vent along the 
highway. The lack of wind causes the drift of the plume from the vent to be very 
unpredictable. It is not unreasonable to expect a decline in the density of several orders 
of magnitude over the distances spanned by the data set. 

Appendix A displays several data collections by AHI over the same road that 
FLYSPEC and COSPEC detected the SO 2 plume. There is no apparent SO 2 in the AHI 
data. There could be several explanations for this. The time of the AHI data collection 
over the road was slightly different than FLYSPEC and COSPEC. The unpredictable 
plume drift due to lack of wind makes it difficult to come to a definite conclusion as to 
whether the plume was present on the road at the time of the AHI data collection. 
Another explanation may be found in the sensitivity of the sensor. As seen in Figure 59 
of Chapter V, there appears to be a lower limit of detection for a down-looking IR sensor 
of 100 to 200 ppm. These values are multiplied by the model’s plume thickness to 
generate concentration path lengths of 15,000 to 30,000 ppm-m. These values are much 
greater than any seen by FLYSEPEC or COSPEC along the road. 

Overall, AHI successfully demonstrated the ability to detect and quantify a 
volcanic SO 2 plume. It is believed that if the experiment would be rerun on a day with 
more ideal weather conditions that may produce data sets that are better geographically 
correlated, a strong comparison of the results could be produced using the methods 
described above. With additional data collections and comparisons with COSPEC and 
FLYSPEC results, AHI could prove to be a valuable addition to the sensors that are 
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currently used to monitor SO 2 plumes. The airborne aspect allows collection of data 
during times that may not be safe for COSPEC or FLYSPEC to collect. It also allows a 
visualization of a large area of the plume nearly instantaneously. With multiple passes 
over a plume, this is could prove beneficial in monitoring the behavior of the plume 
under different weather conditions. The infrared nature of AHI provides the additional 
benefit of monitoring the thermal aspects of the volcano such as active lava flows and 
determining various emissivities of the basaltic rock. Finally, the hyperspectral nature of 
AHI must not be overlooked. A detailed classification of materials, whether it is 
vegetation, rocks, or an SO 2 plume cannot be accomplished without this. 
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APPENDIX A. FALSE COLOR SIMULATED TIMS IMAGES 
WITH DECORRELATION STRETCH APPLIED 



Radiance Image False Color Decorrelation Stretch Image 


Figure 70. Decorrelation Stretch Applied to Simulated TIMS Image 

of R1 A 093134 
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Radiance Image 


False Color Decorrelation Stretch Image 


Figure 71. 


Decorrelation Stretch Applied to Simulated TIMS Image 
of R1A 093226 
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Radiance Image False Color Decorrelation Stretch Image 


Figure 72. Decorrelation Stretch Applied to Simulated TIMS Image 

of R1A 093433 
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Radiance Image False Color Decorrelation Stretch Image 


Figure 73. Decorrelation Stretch Applied to Simulated TIMS Image 

of R1_A_093525 
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Radiance Image False Color Decorrelation Stretch Image 


Figure 74. Decorrelation Stretch Applied to Simulated TIMS Image 

of R1 A 093618 
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Radiance Image 
(Image only 800 lines) 



False Color Decorrelation Stretch Image 
(Image only 800 lines) 


Figure 75. Decorrelation Stretch Applied to Simulated TIMS Image 

of R1 A 093718 
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Figure 76. Decorrelation Stretch Applied to Simulated TIMS Image 

of R1 B 094450 


99 












Radiance Image False Color Decorrelation Stretch Image 


Figure 77. Decorrelation Stretch Applied to Simulated TIMS Image 

of R1 B 094542 
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Radiance Image False Color Decorrelation Stretch Image 

Figure 78. Decorrelation Stretch Applied to Simulated TIMS Image 

of R1 B 094635 
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Radiance Image False Color Decorrelation Stretch Image 


Figure 79. Decorrelation Stretch Applied to Simulated TIMS Image 

of R1_C_095258 
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Radiance Image 


False Color Decorrelation Stretch Image 


Figure 80. 


Decorrelation Stretch Applied to Simulated TIMS Image 
of R1C 095349 
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Radiance Image False Color Decorrelation Stretch Image 


Figure 81. Decorrelation Stretch Applied to Simulated TIMS Image 

of R2 B 100030 
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Figure 82. Decorrelation Stretch Applied to Simulated TIMS Image 

of R2 B 100121 
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Radiance Image 



False Color Decorrelation Stretch Image 


Figure 83. Decorrelation Stretch Applied to Simulated TIMS Image 

of R2 B 100213 
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Radiance Image 


False Color Decorrelation Stretch Image 


Figure 84. 


Decorrelation Stretch Applied to Simulated TIMS Image 
of R2 B 100312 
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Radiance Image False Color Decorrelation Stretch Image 


Figure 85. Decorrelation Stretch Applied to Simulated TIMS Image 

ofR2 B 100403 
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Radiance Image 


False Color Decorrelation Stretch Image 


Figure 86. 


Decorrelation Stretch Applied to Simulated TIMS Image 
of R2 B 100455 
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Radiance Image False Color Decorrelation Stretch Image 


Figure 87. Decorrelation Stretch Applied to Simulated TIMS Image 

of R2 B 100549 
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APPENDIX B. MODTRAN CODE 



Figure 88. MODTRAN Code Used for AHI Data Analysis 
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APPENDIX C 


Figure 89. 


MODTRAN AND AHI RADIANCE PROFILE 
COMPARISONS 


RED CLASS S02: T = 355 K S02 = 475 ppm X:52 Y:635 [5x5 AVG] 



RED CLASS S02: T = 352 K S02 = 350 ppm X:87 Y:534 [5x5 AVG] 




MODTRAN and AHI Radiance Profile Comparison for 
Pixels of the Red SO 2 Class 
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Figure 90. 


RED CLASS S02: T = 350 K S02 = 325 ppm X:88 Y:638 [5x5 AVG] 



RED CLASS S02: T = 361 K S02 = 450 ppm X:115 Y:534 [5x5 AVG] 




Wavelength (p) 


MODTRAN and AHI Radiance Profile Comparison for 
Pixels of the Red SO 2 Class Continued 
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Figure 91. 


RED CLASS S02: T = 346 K S02 = 650 ppm X:39 Y:528 [5x5 AVG] 



RED CLASS S02: T = 335 K S02 = 550 ppm X:154 Y:637 [5x5 AVG] 



RED CLASS S02: T = 345 K S02 = 725 ppm X:51 Y:660 [5x5 AVG] 



MODTRAN and AHI Radiance Profile Comparison 
for Pixels of the Red SO 2 Class End 
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ORANGE CLASS S02: T = 351 K S02 = 625 ppm X:103 Y:646 [5x5 AVG] 



ORANGE CLASS S02: T= 349 K SO2 = 850 ppm X:81 Y:646 [5x5 AVG] 



ORANGE CLASS S02: T = 350 K S02 = 950 ppm X:35 Y:646 [5x5 AVG] 



Figure 92. MODTRAN and AHI Radiance Profile Comparison 
for Pixels of the Orange SO 2 Class 
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ORANGE CLASS S02: T = 343 K S02 = 1400 ppm X:98 Y:513 [5x5 AVG] 



ORANGE CLASS S02: T = 343 K S02= 975 ppm X:66 Y:648 [5x5 AVG] 




Figure 93. MODTRAN and AHI Radiance Profile Comparison 
for Pixels of the Orange SO 2 Class Continued 
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ORANGE CLASS S02: T = 343 K S02 = 1225 ppm X:141 Y:513 [5x5 AVG] 




GREEN CLASS S02: T = 392 K S02 = 575 ppm X:83 Y:522 [4x1 AVG] 



Figure 94. MODTRAN and AHI Radiance Profile Comparison 
for Pixels of the Orange SO 2 Class End and for Pixels of the Green SO 2 Class 


118 




















ORCHID CLASS S02: T = 344 K S02 = 1100 ppm X:140 Y:664 [6x3 AVG] 



ORCHID CLASS S02: T = 323 K S02 = 950 ppm X:210 Y:687 [6x3 AVG] 



ORCHID CLASS S02: T = 317 K S02 = 1275 ppm X:84 Y:671 [6x3 AVG] 



Figure 95. MODTRAN and AHI Radiance Profile Comparison 
for Pixels of the Orchid SO 2 Class 
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CYAN CLASS S02: T = 355 K S02 = 1575 ppm X:38 Y:517 [6x3 AVG] 



CYAN CLASS SQ2: T = 362 K SQ2 = 1600 ppm X:145 Y:511 [6x3 AVG] 



CYAN CLASS S02: T = 371 K S02 = 1150 ppm X:137 Y:521 [6x3 AVG] 



Figure 96. MODTRAN and AHI Radiance Profile Comparison 
for Pixels of the Orchid SO 2 Class 
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BLUE CLASS S02: T = 311 K S02 = 725 K X:199 Y:649 [2x1 AVG] 



BLUE CLASS S02: T = 311 K S02 = 800 ppm X:184 Y:653 [2x1 AVG] 



BLUE CLASS SQ2: T = 316 K SQ2 = 800 ppm X:99 Y:660 [1x1 AVG] 



Figure 97. MODTRAN and AHI Radiance Profile Comparison 
for Pixels of the Orchid SO 2 Class 
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